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Abstract 

We study the chromatic polynomial Pg{q) for mxn triangular-lattice strips 
of widths m < 12p,9f (with periodic or free transverse boundary conditions, 
respectively) and arbitrary lengths n (with free longitudinal boundary condi- 
tions). The chromatic polynomial gives the zero-temperature limit of the par- 
tition function for the g-state Potts antiferromagnet. We compute the transfer 
matrix for such strips in the Fortuin-Kasteleyn representation and obtain the 
corresponding accumulation sets of chromatic zeros in the complex g-plane in 
the limit n ^ oo. We recompute the limiting curve obtained by Baxter in the 
thermodynamic limit m, n — > oo and find new interesting features with possi- 
ble physical consequences. Finally, we analyze the isolated limiting points and 
their relation with the Beraha numbers. 



Key Words: Chromatic polynomial; chromatic root; antiferromagnetic Potts model; 
triangular lattice; transfer matrix; Fortuin-Kasteleyn representation; Beraha-Kahane- 
Weiss theorem; Beraha numbers. 



1 Introduction 



The antiferromagnetic g-state Potts model [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] 
exhibits unusual behavior not found in ferromagnets. Indeed, its phase diagram and 
critical behavior depend crucially on the lattice structure, unlike the ferromagnetic 
case where the concept of universality applies. Specifically, for each lattice C there 
exists a number qc{C) such that for all q > Qc the model is disordered at any temper- 
ature, including zero temperature [11]. Exactly at g = Qc{j^) the system is disordered 
at all positive temperatures and has a zero-temperature critical point. The zero- 
temperature limit of the antiferromagnetic Potts model is particularly interesting 
because its partition function on a finite graph G coincides with the chromatic poly- 
nomial Pg{<1)i which counts the number of ways of coloring the vertices of G using q 
colors subject to the constraint that adjacent vertices always receive different colors 
[15]. 

In Refs. [16, 17] we undertook a study of the zeros of the chromatic polynomial 
when the parameter q is allowed to take complex values (see [16] for detailed references 
to the previous literature). In those papers we studied strips of the square lattice; 
here we extend that work to the triangular lattice. The triangular-lattice case is 
of particular interest because the path-breaking work of Baxter [18, 19] provides a 
conjectured exact solution in the thermodynamic limit. 

The study of the complex zeros of the chromatic polynomial is inspired by the 
Yang-Lee picture of phase transitions [20] . We study families of graphs Gn for which 
the chromatic polynomial can be expressed via a transfer matrix of fixed size M x M: 

PgM) = iM{q)Tm (1-la) 

M 

= 5]«fe(g)Afc(gr , (1.1b) 

fc=i 

where the transfer matrix T{q) and the boundary- condition matrix A{q) are polyno- 
mials in g, so that the eigenvalues {A/,.} of T and the amplitudes {a^} are algebraic 
functions of q. Rather than using T(g) to compute the zeros of the chromatic poly- 
nomial for a finite strip m x n, we have focussed on the direct calculation of their 
accumulation points in the limit n — >■ oo, i.e. for the case of an semi-infinite strip 
[21, 22, 23, 24, 16, 17]. According to the Beraha-Kahane- Weiss theorem [25, 26, 27], 
the accumulation points of zeros when n ^ oo can either be isolated limiting points 
(when the amplitude associated to the dominant eigenvalue vanishes, or when all 
eigenvalues vanish simultaneously) or belong to a limiting curve B (when two dom- 
inant eigenvalues cross in modulus). By studying the limiting curves for different 
values of the strip width m, we hope to learn new features of the thermodynamic 
limit m — s> oo. 

To determine the isolated limiting points, we shall take advantage of the following 
simple result [22]: 

M 

det D = l[ak n (-^^-^0', (1-2) 

k=l l<i<j<M 
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where D is the M x M matrix with entries Dij = X]fc=i 0(k{^ky~^-'~'^ — PGi+j-2- 

An important feature of the hmiting curve B is that it typically crosses the positive 
real axis at a point qo{m)} Physically, qoim) corresponds to a point of non-analyticity 
of the ground-state degeneracy per site [23, 28]. As the strip width m grows, this 
crossing point qo{m) increases and presumably tends to a limiting value 50(00) ■ On 
the other hand, as m ^ cxd the curve B = Bm presumably tends to a thermodynamic- 
limit curve Boo- We define to be the largest value where B^o crosses the real g-axis: 
please note that qc may or may not correspond to crossings of the real axis for any 
finite m. In general the value qo{oo) is smaller than q^., although for some lattices 
they may coincide (this depends on the shape of the curves Bm and Boo). Indeed, in 
Ref. [17] evidence was presented suggesting that for the square lattice go (00) = gc = 3. 

The crucial role of go is further emphasized by studying the relation between 
chromatic polynomials and the so-called Beraha numbers 

S„ = 4cos^- = 2 + 2COS— forn = 2,3, ... . (1.3) 
n n 

It has been foTind in a number of cases [29, 21, 22] that chromatic roots tend to 
accumulate at some of the Beraha numbers. In the case of the square lattice, we 
have found empirically [16, 17] that on a strip of width m with either free or periodic 
transverse boundary conditions^, there is at least one vanishing amplitude ai{q) at 
each of the first m Beraha numbers S2, . . . , S^+i (but not higher ones). It thus 
appears that in the limit m — > 00 all the Beraha numbers will be zeros of some 
amplitude.'^ Moreover, we found that the vanishing amplitude corresponds to the 
eigenvalue obtained by analytic continuation in g from the one that is dominant 
at small real g, in agreement with a conjecture of Baxter [19, p. 5255]. Thus, the 
first few Beraha numbers — namely, those (up to at most S^+i) that lie below 
the point qo{m) — correspond to the vanishing of a dominant amplitude and hence, 
by the Beraha-Kahane-Weiss theorem, to a limit point of chromatic roots, while 
the remaining Beraha numbers do not. As the strip width m grows, qo{m) tends 
to go (00), and the limiting points of chromatic roots are thus constrained to be the 
points B2, B2, . . . ,Bp < go(oo). This scenario for the accumulation of chromatic roots 
at some of the Beraha numbers was set forth by Baxter [19] and elaborated by Saleur 
[8]; further references can be found in [16]. 

^If there is more than one such crossing, we define qoim) to be the smallest such crossing. When 

no such crossing occurs, the Hmiting curve often inchidcs a pair of complex-conjugate cndpoints 
rather close to the positive real q-axis. In these cases, we define qo{m) to be the point closest to 
that axis with positive imaginary part. 

^Let m (resp. n) denote the number of sites in the transverse (resp. longitudinal) direction of 

the strip, and let F (resp. P) denote free (resp. periodic) boimdary conditions in a given direction. 
Then we use the terminology: free (mp x np), cylindrical (mp x np), cyclic (top x np), and toroidal 
(top X np). In this paper we consider free and cylindrical boundary conditions, as well as a new 
type of boundary condition that we shall call "zig-zag" (Section 5). 

^As we shall sec, a similar (but not identical) statement appears to hold true also for the triangular 
lattice: see Section 7.2. 
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In the present publication we shall be concerned with the antiferromagnetic Potts 
model on the triangular lattice. For this case, Baxter and collaborators [30, 18, 19] 
have determined the exact free energy (among other quantities) on two special curves 
in the {q, -plane: 

v^ + 3v^-q = (1.4) 
V ^ -1 (1.5) 

The uppermost branch [v > 0) of curve (1.4) is known to correspond to the ferro- 
magnetic critical point [30, 3] , and Baxter [18] initially conjectured (following a hint 
of Nienhuis [31]) that the zero-tempcraturc antiferromagnetic model (1.5) is critical 
in the interval < g < 4. This prediction is known to be correct for q = 2 [32, 33, 34] 
and is believed to be correct also for g = 4 [5, 35, 36]. On the other hand, for g = 3 
the conjecture contradicts the rigorous result [37], based on Pirogov-Sinai theory, that 
there is a low-temperature phase with long-range order and small correlation length.^ 
In any case, for g > 4 we expect that the triangular-lattice Potts model is noncritical 
even at zero temperature; this has recently been confirmed by Monte Carlo simulation 
of the models with g = 5, 6 [38]. We therefore expect that for the triangular lattice 
gc = 4. 

For the model (1.5), Baxter [18] used a Bethe Ansatz to compute three different 
expressions gi{q) [i = 1, 2, 3] that he argued correspond to the dominant eigenvalues of 
the transfer matrix in different regions Vi of the complex g-planc; in a second paper 
[19] he provided corrected estimates for the precise locations of T>i,T>2,'D^. Using 
these formulae, he determined the value of go(oo) as 

go(oo,tri) 3.81967 (Baxter) (1.6) 

One important outcome of the present paper (see Section 6.2) is that the "phase 
diagram" predicted by Baxter's formulae is actually more complicated than what 
Baxter found. In particular, it now appears that the correct value of go(cxo, tri) coming 
from Baxter's eigenvalues Qi is slightly smaller than Baxter's value (1.6), the corrected 
value being 

go(oo, tri) = S12 = 2 + 3.73205 (this paper) (1.7) 

We therefore conjecture that the isolated limiting chromatic roots of the infinite- 
size triangular lattice are B2, . . . , Bu and possibly B12, rather than B2, . . . , Bn as 
conjectured by Baxter. 

To study the approach to the thermodynamic limit, we have computed the transfer 
matrix for triangular-lattice strips of widths 2 < m < 9 with free boundary conditions 
and 2 < m < 12 with cylindrical boundary conditions, and we have determined the 
corresponding limiting curves B. We have also undertaken a detailed comparison of 
our finite-lattice eigenvalues with Baxter's eigenvalues Qi, and have found a surprising 

"'A Monte Carlo study of the g = 3 model found strong evidence for a first-order transition to an 
ordered phase at /3J w -1.594 [10]. 
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fact (Section 6.3): in the region of the g-plane corresponding to Baxter's V^, we find 
no evidence in our finite-lattice data of any eigenvalue corresponding to gi. As a 
consequence, the "extra" branches of the phase diagram found in Section 6.2 appear 
to be absent after all! There seem to be two possibihties: 

(a) The eigenvalue gi really is present in the region P3, but only for strip widths 
much larger than those we have studied. In this case, the limiting curve Boo 
really would exhibit all the complexities found in Section 6.2, and the correct 
value of qo would be given by (1.7) rather than (1.6). 

(b) For some reason, the eigenvalue gi is not present in this region (though it is 
clearly present elsewhere). In this case, the limiting curve Boo would be similar 
to that depicted by Baxter, and the correct value of qo would be given by (1.6) 
after all. 

We discuss these issues further in Section 6.4. 

Our finite-lattice data also serve as a testing ground for the general conjectures 
on Beraha numbers as stated above (see Ref. [16] for further details). We discuss this 
further in Section 7.2. 

Previous studies using a similar transfer-matrix approach have been made notably 
by Shrock and collaborators. In particular, they have considered triangular-lattice 
strips of width m < 5 for free and cylindrical boundary conditions [24, 41, 42, 43, 44].^ 
They have also considered other boundary conditions for the same lattice [44, 45, 
46, 47].^ Generalizations to nonzero temperature for several boundary conditions 
have been carried out in Refs. [48, 49]. Finally, Refs. [47, 50] discuss some general 
structural properties of the Potts-model partition function and chromatic polynomial 
on square-lattice and triangular-lattice strips. 

This paper is laid out as follows: In Section 2 we discuss some brief preliminaries. 
In Section 3 we give our numerical results for free transverse boundary conditions, and 
in Section 4 for periodic transverse (cylindrical) boundary conditions. We have also 
examined a third type of boundary conditions, called "zig-zag" , which we introduce 
and motivate in Section 5. In Section 6 we analyze Baxter's [18, 19] exact solution 
for the thermodynamic limit and carefully recompute his phase diagram, finding 
interesting new features with possible physical consequences; we also compare his 
predictions for the dominant eigenvalues with our finite-lattice data and comment on 
the agreements and discrepancies. Finally, in Section 7 we present our conclusions. 

^The case m = 4 with cylindrical boundary conditions was first done by Beraha and Kahane [21]. 
The case m = 5 with cylindrical boundary conditions was first done by Beraha, Kahane and Weiss 
[22]. 

^The case 2p x np (i.e., cyclic boundary conditions) was first done by Beraha, Kahane and Weiss 
[22]. 
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2 Preliminaries 



The general theory of the transfer-matrix method used here has been explained 
in Ref. [16], and the implementational details of an improved algorithm have been 
given in Rcf. [17]. Suffice it here to mention that we have used the Fortuin-Kasteleyn 
representation [39, 40] of the Potts model in the computation of the transfer matrix; 
therefore, all quantities are expressed as polynomials in q. 

To compute the limiting curves B we have used two different techniques: the 
resultant method and a direct-search method. These techniques have been described 
in [16, Section 4.1], and we use here the same conventions and notation. 

Let us briefly mention a few improvements/additions to our methodology: 

Computation of T points. We have adopted an improved method for locating T 
points, based on applying a numerical minimization algorithm (e.g. Mathematica's 
FindMinimum) to the function 

F{q) = - \H<i)\f + (|Ai(g)| - \\M\f + (|A2(?)| - |A3(g)|)^ , (2.1) 

where Ai, A2, A3 arc the three eigenvalues of the transfer matrix of largest modulus. 
At any given q these eigenvalues can easily be computed numerically by finding the 
roots of the characteristic polynomial of the transfer matrix. Using this method, we 
are able to locate T points far more precisely than in our previous work. 

Fixed zeros. When q is an integer and the graph G is not g-colorable, we have 
Pg{<1) = 0. For this reason, certain small integers q can be "fixed" zeros of the 
zero-temperature partition function, independent of the strip length n. In particular, 
g = 0, 1 are roots for all widths m > 2 and lengths n>l. Furthermore, = 2 is a 
root for all triangular-lattice strips of widths m > 2 and lengths n > 2, because the 
triangular lattice is not bipartite. Finally, g = 3 is a root for all cylindrical triangular- 
lattice strips of widths that are not multiples of 3 (with lengths n > 2), because these 
graphs are not 3-colorable. 

It is interesting to see how these behaviors come about from the point of view of 
the transfer-matrix formalism. The partition function on a lattice of length n has the 
form 

M 

Zn = $^«,(g)A,(g)'^-\ (2.2) 

where the {Xk\ are the eigenvalues of the transfer matrix and the {0;^} are the cor- 
responding amplitudes. A particular value q can then be a "fixed" zero of Z for any 
of three reasons: 

1) All the amplitudes ctfe vanish at q. Then Zn{q) = for all n > 1. 

2) All the eigenvalues A^ vanish at q. Then Z„(g) = for all n > 2. 

3) "Mixed case": Neither all the amplitudes nor all the eigenvalues vanish at g, 
but for each k either ak or A^ vanishes at q (or both). Then Zn{q) = for all 
n>2. 
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As we shall see, the points g = and q = 1 will be fixed roots belonging to Case 1: 
all the amplitudes vanish due to an overall prefactor q{q — 1). The point q = 2 will 
be a fixed root belonging sometimes to Case 1 and sometimes to Case 3 (and to Case 
2 when m = 2 for all boundary conditions) . For cylindrical strips where the width is 
not a multiple of 3, the point q — 3 will be a fixed root belonging to Case 2 for m = 4p 
and to Case 3 for m > 5p. We shall endeavor to explain in each case the mechanism 
underlying the fixed zeros; these results will be summarized in Section 7.3. 

Computation of isolated limiting points. To find the isolated limiting points, we 
first compute symbolically the determinant det D{q) defined by (1.2); this determinant 
is a polynomial in q with integer coefficients, typically of very high degree. We 
then compute numerically the zeros of this polynomial, using the MPSolve 2.1.1 
package [51, 52]; these zeros correspond to points q where at least one amphtude 
ak{q) vanishes. Finally, we test numerically each of these zeros to see whether the 
amplitude corresponding to the dominant eigenvalue is vanishing; if it is, then the 
point in question is an isolated limiting point. This method is guaranteed to discover 
all of the isolated hmiting points. We shall not bother to report here all the zeros 
of det D{q), but only (a) the isolated hmiting points and (b) the Beraha numbers S„ 
that are zeros of some subdominant amplitude. Please note that whenever Bn is a 
zero of det D{q), so are all the primitive generalized Beraha numbers 



where k is relatively prime to n, since they have the same minimal polynomial Pniq) 
[16, Section 2.3]. 

Unfortunately, in some cases the matrix D{q) is so large that we have been unable 
to compute symbolically its determinant. In these cases, it is more convenient to 
compute numerically the eigenvalues {Xj} and their corresponding amplitudes {aj} 
and check (a) whether any of the amplitudes vanish and (b) whether the amplitude 
a* associated to the dominant eigenvalue vanishes. We have restricted our search to 
certain "candidate" values of q (or neighborhoods in the complex g-plane), namely 
(a) the Beraha numbers i?„ for n < 50, and (b) any real or complex values oi q ^ B 
where zeros of seem to be accumulating as n gets large. When there is a exact 
candidate (such as the Beraha numbers S„), we have computed the amplitudes with 
high-precision arithmetic (200 digits of precision at least). We considered that an 
amplitude is zero when its absolute value is less than (for instance) 10^^^°. When we 
do not have an exact candidate, we tried to minimize the dominant amplitude \a^\ 
around a region where zeros of Z„ tend to accumulate as n grows. This situation 
occurred only for the strips of widths 8z and lOz ("zig-zag" boundary conditions, see 
Section 5). The condition |q;*| ^ 10~^^ holds for all the cases reported here. 

In all the cases where we are unable to compute det D{q) symbolically, we are able 
to assert that certain points are indeed isolated limiting points, but we cannot claim 
with confidence that we have found all of the isolated limiting points. 





2 + 2 cos 



n 



(2.3) 
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3 Numerical Results for the Triangular-Lattice 
Chromatic Polynomial: 
Free Boundary Conditions 

We have computed the transfer matrix T{mp) and the hmiting curves B for 
triangular-lattice strips of widths 2 < m < 9 with free boundary conditions in both 
directions. We also write as a synonym for the strip width m. 

As explained in Ref. [16], the chromatic polynomial for this family of strip lattices 
can be written as 

Z{mF X np) = u^HT(mF)""Vid (3.1) 

where u and Vid are certain vectors, and T^mp) = VH is the transfer matrix. Here 
H (resp. V) corresponds to adding one new layer of horizontal (resp. vertical and 
diagonal) bonds: sec Figure 1(a). The matrices H, V and T{m-F) act on the space of 
connectivities of sites in the top layer, whose basis elements v-p are indexed by parti- 
tions V of the single-layer vertex set {1,2,..., m}. In particular, Vid = V{{i}^{2},...,{m}}- 
Since the strip lattices we are dealing with are planar, only non-crossing partitions 
V can occur. 

In the particular case of the chromatic polynomial (i.e. the zero-temperature an- 
tiferromagnet), the horizontal operator H is a projection (H^ = H), and we can work 
in its image subspace by using the modified transfer matrix T'{mp) — HVH in place 
of T(mF) = VH, and using the basis vectors 

w-p = Hvp (3-2) 

in place of v-p. Then we can rewrite (3.1) as 

Z{mF X np) = u^T'(mF)""^Wid (3.3) 

where Wjd = Hvjd. Please note that wp = Hvp = for any partition V that includes 
nearest-neighbor sites in the same block, so we can ignore all such partitions. The 
dimension of the transfer matrix T'(mF) is therefore equal to the number of non- 
crossing non- nearest- neighbor partitions of the set {1,2, .. . ,m}, which is given by 
the Motzkin number Mm-i [16]. To simplify the notation, we will drop the prime 
in T'(mF) and denote the basis elements wp by a shorthand Tising Kronecker delta 
functions: for instance, W|{i 3} {2},{4,6},{5}} will be written di^d^Q. We denote the set 
of basis elements for a given strip as P = {wp}. For instance, the basis for m = 3 is 

P = {l,5l3}- 

We have checked the self-consistency of our finite-lattice results using the trivial 
identity 

Z{m-p X np) ~ Z{nF x mp) (3.4) 

for all pairs 2 < m, n < 9. 
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3.1 Lx — 

This case is trivial, as the transfer matrix is one-dimensional: 

Z(2f xnp) = g(g-l)(g-2)2("-i) (3.5) 

Since there is only one eigenvalue, there is obviously no crossing, hence B — 0. 
However, there are zeros for all n at g = 0, 1 and for n > 2 at g = 2. The fixed zeros 
at g = 0, 1 arise from a vanishing amplitude, and the fixed zero at g = 2 arises from 
a vanishing eigenvalue. 

3.2 Lx — 3p 

The transfer matrix is two-dimensional. In the basis P = {1, Sis} it can be written 

as 

^ / g^-7g2 + 17g-14 g2-6g + 9 \ 
and the partition function is equal to 

T 

Z(3p X np) = g(g - 1) ^ 7 ^ ) • r(3p)--i • J ) (3.7) 

The limiting curve B (sec Figure 2) contains three disconnected pieces and it 
crosses the real axis at go ~ 2.5698402910. There are six endpoints: 

g 1.2047381150 ± 1.1596169599 i (3.8a) 
g ^ 2.3930361082 ± 0.2538745688 i (3.8b) 
g ^ 3.4022257768 ± 0.5865084714 i (3.8c) 

These results were previously obtained by Rocek et al. [24]. 
The determinant detL'(g) has the form 

deiD{q)^q\q-l)\q-2f (3.9) 

Thus, it vanishes at the first three Beraha numbers g = 0, 1, 2. At those points the 
dominant amplitude vanishes, hence they are isolated limiting points. 

In fact, the partition function vanishes at g = 0, 1 for all n, and at g = 2 for all 
n > 2. Obviously, at g = 0, 1 both amplitudes vanish, due to the prcfactor g(g — 1) 
in (3.7). [This happens for all strips of width m > 2; we will henceforth call these 
zeros "trivial".] At the fixed zero g = 2, there is one nonzero eigenvalue (A* = — 1) 
with a vanishing amplitude and one zero cigcnvahic with a nonvanishing amplitude; 
we are therefore in Case 3 described in Section 2. The fourth real zero (see Table 1) 
converges at an approximate 1/n rate to the value go ~ 2.5698402910. 
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3.3 Lx = 4f 



The transfer matrix is four-dimensional. In the basis P = {1, ^13, ^24, (^14}, it takes 
the form 



T(4p) 



Til 

-q^ + bq 
-q^ + 5q 
q-2 



T12 

— 6 — 5g + 6 

— 6 — 5g + 6 



Ti3 Tu \ 

-q + 3 -2(g-3) 

-6g + 9 g2_8g+i5 



(3.10) 



g^ - 5g + 6 -g + 3 g^ - 7g + 13 / 



where 



Til 
T12 

Tl3 
Ti4 



Q - 



Q - 



lOg^ + 39g2 - 70g + 48 
9g2 + 26g - 24 
9g2 + 28g - 30 
lOg^ + 36g - 45 



(3.11a) 
(3.11b) 
(3.11c) 
(3.11d) 



The partition function is equal to 



Z(4f X np) = g(g - 1^ 



g-1 
g-1 
g-2 



■T(4j 



/ 1 \ 






(3.12) 



The limiting curve B (see Figure 3) contains two complex-conjugate disconnected 
pieces that do not cross the real axis. The closest points to the real axis are go ~ 
2.7592502040 ± 0.1544431251 i. There are six endpoints: 



Q 



0.8164709452 ± 1.2804094073 i 
2.7592502040 ± 0.1544431251 i 
3.6398304896 ± 0.5986827987 i 



(3.13a) 
(3.13b) 
(3.13c) 



There are T points at g fa 3.3341785562 ± 0.8829730283 i. These results were previ- 
ously obtained by Rocek et al. [24, 41]. 
The determinant det -D(g) has the form 

detL>(g) = -g^(g-l)^(g-2)i3(g2-3g+l)(g-3)^(g^-llg^ + 46g2-86g+61)2 (3.14) 

Wc recognize, as factors in (3.14), the first five minimal polynomials Pk{q) for the 
Beraha numbers Bk [16, Table 1]; therefore dctD(g) vanishes at the first five Bcraha 
numbers g = 0,l,2,i?5,3. The dominant amplitude vanishes at g = 0, 1, 2, i?5, so 
these are isolated limiting points. At g = 3 only two subdominant amplitudes vanish, 
so this is not an isolated limiting point. Similarly, all of the zeros of the last factor in 
(3.14) correspond to the vanishing of subdominant amplitudes only, so none of them 
is an isolated limiting point. 

In fact, the partition function vanishes at g = 0, 1 for all n, and at g = 2 for 
all n > 2. The fixed zeros at g = 0, 1 are trivial. At g = 2, there are three zero 
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eigenvalues and a unique leading eigenvalue A* = 4 with zero amplitude. Notice that 
the transfer matrix is not diagonal for q = 2: there is a 2 x 2 nontrivial Jordan block 
corresponding to A = and whose contribution to the partition function is zero for all 
n. The amplitude corresponding to the other A = eigenvalue is 2; we are therefore 
in Case 3 described in Section 2. Finally, the fourth real zero converges exponentially 
fast to i?5 (see Table 1). 

Please note that for this strip there is a vanishing subdominant amplitude at 
q — Bk for a Beraha number k greater than m + 1 (namely, Bq = 3). As we shall 
see, this occurs frequently for the triangular lattice, and contrasts with the behavior 
observed for the square lattice [16, 17]. 

3.4 Lx = 5f 

The transfer matrix is nine-dimensional; it can be found in the Mathematica file 
transfers, m available as part of the electronic version of this paper in the cond-mat 

archive. This strip has been previously studied by Chang and Shrock [44]; but they 
did not compute the limiting curve. 

The limiting curve B is connected (see Figure 4). It crosses the real axis at qo = 3. 
There are six endpoints: 

q ^ 0.5586170364 ± 1.2816149610 i (3.15a) 
q ft! 3.0474871745 ± 0.8171660680 i (3.15b) 
q ft 3.7782975917 ± 0.5699779858 i (3.15c) 

The topology of the hmiting curve is rather involved. It has 12 T points: q ft 
3.1572589261±0.7951215102i, q ft 3.1251751109±0.8152460413 i, q ft 3.2093444343± 
0.9296294663 q ft 3.3452062643±0.9758086833i, q ft 3.3248793469±0.9987588766 
and q ft 3.2492362818 ± 1.11850738091 These points define four closed regions. The 
first five T points are the vertices of two complex-conjugate regions which look ap- 
proximately like rectangular bands. The third, fifth and sixth T points above define 
two complex-conjugate triangular-like regions. 

The determinant det D{q) is given by 

det D{q) = q%q - ifiq - 2f\q^ - 3g + l)\q - 3f'P{qf (3.16) 

where the polynomial P{q) can be found in the Mathematica file transf erS .m. 
The factors appearing in det D{q) are the first five polynomials Pk{q) given in [16, 
Table 1]; therefore det D{q) vanishes at g = 0, l,2,i?5,3. The dominant amplitude 
vanishes at the first four of them (but not at g = 3), so that g = 0, 1, 2, are isolated 
limiting points. All of the zeros of P{q) correspond to the vanishing of subdominant 
amplitudes only, so none of them is an isolated limiting point. 

The first two real zeros g = 0, 1 are trivial ones. At g = 2 there are two nonzero 
eigenvalues with zero amplitudes and seven zero eigenvalues. The transfer matrix is 
not diagonalizable: there are two Jordan blocks of dimension 3 and 2, respectively, 
associated to the eigenvalue A = 0. The contribution of these blocks to the partition 
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function vanishes for all n. In addition, the amplitude corresponding to the other two 
zero eigenvalues are and 2; we are therefore in Case 3 described in Section 2. The 
fourth real zero converges exponentially fast to the value q = (see Table 1); and 
the fifth real zero converges at an approximate 1/n rate to the value go = 3- This 
agrees with the fact that g = 3 is a regular limiting point and not an isolated limiting 
point. 

3.5 Lx = 6f 

The transfer matrix is 21-dimensional; it can be found in the Mathematica file 
transf er3 .m. 

The limiting curve B is connected (see Figure 5). It crosses the real axis at 
qo ~ 3.1609256737. There are six endpoints: 

q ^ 0.3796307748 ± 1.2450702104 i (3.17a) 
q fti 2.9641235697 ± 1.1179839989i (3.17b) 
q ^ 3.8664092416 ± 0.5329463088 i (3.17c) 

There are four T points at g 3.3081144403±1.2171494282 i, and q 3.5005856709± 
0.9442298756 i. 

The determinant det D{q) is given by 

detD{q) = q^\q-lf\q-2f''\q^-3q + iy%q-3)'''%q^-5q^ + 6q-l)P{q)^ (3.18) 

where the polynomial P{q) can be found in the Mathematica file transf erS.m. 
The factors appearing in det D{q) are the first six polynomials Pk{q) given in [16, 
Table 1]. Thus, det D{q) vanishes ai q — 0, 1, 2, S5, 3, i^y. The dominant amplitude 
vanishes at the first five of them (but not at Sy), so that g = 0, 1, 2, B5, 3 are isolated 
limiting points. All of the zeros of P{q) correspond to the vanishing of subdominant 
amplitudes only, so none of them is an isolated limiting point. 

The first two real zeros q = 0,1 are trivial ones. At g = 2 we get six nonzero 
eigenvalues with zero amplitudes and 15 zero eigenvalues. We find again that the 
transfer matrix is not diagonalizable for g = 2: there are five nontrivial Jordan blocks 
(one of dimension 3 and four of dimension 2) corresponding to the eigenvalue A = 0, 
and whose contribution to the partition function always vanishes. The amplitudes of 
the other four zero eigenvalues are (0,0,0,2). This combination seems to be the generic 
case for free boundary conditions: all amplitudes vanish except one corresponding to 
a zero eigenvalue. When the transfer matrix is not diagonalizable, then the nontrivial 
Jordan blocks correspond to A = and their contribution is always zero. The fourth 
and fifth real zeros converge exponentially fast to the values q = and g = 3, 
respectively (see Table 1); the sixth real zero converges at an approximate 1/n rate 
to the value go ^ 3.1609256737. 

3.6 Lx = 7f 

The transfer matrix is 51-dimensional; it can be found in the Mathematica 
file transf er3 .m. In this case we have been unable to compute symbolically the 
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resultant, hence the computation of the hmiting curve has been performed using the 
direct-search method. 

The hmiting curve B is connected (see Figure 6). It crosses the real axis at 
go ~ 3.2764013231. There are six endpoints: 

q w 0.250538 ± 1.196864i (3.19a) 
q ^ 3.925804 ± 0.496672 i (3.19b) 
q fa 2.878928 ± 1.343851 i (3.19c) 

There are four T points at q ^ 3.6146786603±0.9081562491 i and g fa 3.2704141478± 
1.5194310419 i. 

We have been unable to compute the determinant det D{q). However, we com- 
puted the amplitudes numerically at each of the Beraha numbers up to B^q and 
determined in particular whether it is an isolated limiting point or not. As always, 
g = 0, 1 are trivial isolated limiting points where all the amplitudes vanish. The 
dominant amplitude also vanishes at g = 2, S5, 3, S7, so they are isolated limiting 
points too. Finally, a subdominant amplitude vanishes at g = Bg. 

The first two real zeros g = 0, 1 arc trivial ones. At g = 2, all amplitudes vanish 
except one corresponding to a zero eigenvalue. In particular, there are 12 nonzero 
eigenvalues with zero amplitudes, 8 zero eigenvalues with zero amplitudes, and one 
zero eigenvalue with nonzero amplitude. We also find 12 nontrivial Jordan blocks 
corresponding to A = 0, whose contribution to the partition function is always zero. 
The fourth, fifth and sixth real zeros converge exponentially fast to the values i?5, 3 
and Bj, respectively (see Table 2); however, the convergence to i?7 ~ 3.246979603717 
is slowed by its nearness to the regular hmiting point go ~ 3.2764013231. For lengths 
n ^ 77, a seventh real zero appears: it converges (at an approximate 1/n rate) to go. 

3.T Lx = 8f 

The transfer matrix is 127-dimensional; it can be found in the Mathematica file 
transfers. m. Again we have used the direct-search method to locate the points of 
the limiting curve B. 

The limiting curve B is connected (see Figure 7). It crosses the real axis at 
go ~ 3.3610599515. There are six endpoints: 

g 0.154432 d= 1.146669 i (3.20a) 
g fa 2.793496 ± 1.521468 i (3.20b) 
g fa 3.967566 ± 0.463648 i (3.20c) 

There are four T points at g 3.2555859898± 1.7000353877 i and g ^ 3.6703287722± 

0.8845072864 i. 

We were unable to compute the determinant dct D{q). However, we computed the 
amplitudes numerically at each of the Beraha numbers Bn up to -B50 and determined 
in particular whether it is an isolated limiting point or not. As always, g = 0, 1 
are trivial isolated limiting points where all the amplitudes vanish. The dominant 
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amplitude vanishes also sd, q — 2, 3^,3, B7, so they are isolated limiting points too. 
Finally, some subdominant amplitudes vanish at q — Bg,Bg; they are not isolated 
limiting points. 

The first two real zeros g = 0, 1 are trivial ones. At the third real zero q — 2, all 
amplitudes vanish except one corresponding to a zero eigenvalue. At g = 2 the transfer 
matrix is not diagonalizable: we find 30 nontrivial Jordan blocks corresponding to 
A = 0, but their contribution to the partition function vanishes for all n > 1. We 
also get 33 nonzero eigenvalues with zero amplitudes, 20 zero eigenvalues with zero 
amplitudes, and one zero eigenvalue with nonzero amplitude. The fourth, fifth and 
sixth real zeros converge exponentially fast to the values B^, 3 and B-^, respectively 
(sec Table 2). We also expect a seventh real zero converging (at an approximate 1/n 
rate) to go ~ 3.3610599515. Such zero docs not appear up to lengths n = 96 (see 
Table 2). We would need to go to very large n to observe this additional zero. 

3.8 Lj; = 9f 

The transfer matrix is 323-dimensional; it can be found in the Mathematica 
file transf er3 .m. The size of this transfer matrix prevented us from computing the 
limiting curve B. However, we were able to compute the point where the limiting 
curve crosses the real axis: it is go ~ 3.4251304673. In Figure 8 we show the zeros of 
Z for the finite lattices Qp X 45? and 9f x 90f. 

The first two real zeros g = 0, 1 arc trivial ones. We have checked numerically 
that the dominant amplitude vanishes at g = 2, B^, 3, Bj, B^ (and at no other Beraha 
numbers up to -B50); these are therefore isolated limiting points. By inspection of 
Figure 8 we do not find any candidate for a non-Beraha real isolated limiting point or 
for a complex isolated limiting point. A subdominant amplitude vanishes for q — Bg 
and BiQ. 

At g = 2 the transfer matrix is not diagonalizable: we find 76 nontrivial Jordan 
blocks (up to dimension 5x5) corresponding to A = 0; none of these blocks contribute 
to the partition function for any n> 1. We also get 74 nonzero eigenvalues with zero 
amplitudes, 50 zero eigenvalues with zero amplitudes, and one zero eigenvalue with 
nonzero amplitude. The fourth, fifth, and sixth real zeros converge exponentially fast 
to the values B^, 3 and By, respectively (see Table 2). We expect that the seventh real 
zero will converge exponentially fast to B^ — 2 + \/2 and that there will be (for large 
enough length n) an additional real zero that will eventually converge to the value 
g ^ 3.4251304673 > Bg at an approximate 1/n rate. However, we would probably 
need to go to very large lengths n in order to observe this behavior. 
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4 Numerical Results for the Triangular-Lattice 
Chromatic Polynomial: 
Cylindrical Boundary Conditions 

We have also computed the transfer matrix T(mp) and the hmiting curves B 
for triangular-lattice strips of widths 2 < m = L^- < 12 with cylindrical boundary 
conditions, i.e. periodic b.c. in the transverse direction and free b.c. in the longitudinal 
direction. 

The partition function can be written analogously to (3.3) as 



where T(mp) = HVH. (See [16, Section 3.1] for some remarks about the actual 
computation of this transfer matrix.) Since T(mp) commutes with translations (due 
to the periodic b.c.) and the vectors u and Wid are translation-invariant, we can 
restrict attention to the translation-invariant subspace. The dimension of T(mp) 
is therefore given by the number of equivalence classes modulo translation of non- 
crossing non- nearest-neighbor partitions. This number is denoted by TriCyl(m) in 
[16, Table 2]. A general analytic expression for TriCyl(m) is not known, although 
such a formula has been obtained for prime values of m [53, Theorem 3]. 

We have checked our results for widths 2 < mp < 12 and lengths rip = 2, 3, 4 by 
comparing to the results of Beraha-Kahane- Weiss [22] (resp. Chang-Shrock [44]) for 
width np = 2 (resp. rip = 3,4) and length mp (i.e. cyclic boundary conditions), using 
the trivial identity 



This is a highly non-trivial check of the correctness of our results. 

The triangular lattice with cylindrical boundary conditions possesses a curious 
reflection symmetry that wc shall now explain. Note first that the triangular lattice is 
not invariant under reflection in the transverse direction^, since reflection changes the 
direction of the diagonal bonds. Nevertheless, in the translation-invariant subspace 
the transfer matrix does commute with reflection, because by translating the upper 
row (of a pair of rows) by one unit, one can change I\ into /[ , thereby converting 
the triangular lattice into a reflected triangular lattice! Because the transfer matrix 
(in the translation-invariant subspace) commutes with reflection, we can pass to a 
new basis consisting of connectivities that are either even or odd under reflection.^ 

'^More precisely, if the strip width is even, one can choose to reflect either through a pair of lattice 
sites or through a pair of bisectors (sites with half- integer coordinates); the two choices differ by 
a subsequent translation. If the strip width is odd, then every reflection axis passes through one 
lattice site and one bisector. 

^Lct us call {fjljli the connectivity basis in the translation-invariant subspace. Some of these 
basis elements are invariant under reflection; the rest can be grouped into pairs {va,Vjs) that map 
into each other under reflection. A basis for the reflection-even (i.e. reflection-invariant) subspace 
is then given by the basis elements in the first set together with the combinations Va + from the 
second set. A basis for the reflection-odd subspace is given by the combinations Va — Vjj from the 
second set. 
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In this new basis, the transfer matrix T(mp) is block-diagonal: 

where T_|_ (resp. T_) corresponds to the reflection-even (resp. reflection-odd) sub- 
space. Furthermore, the reflection-odd components of the start vector Wid = Hvjd are 
identically zero, since both Vid and H are manifestly reflection- invariant. Likewise, 
the reflection-odd components of the final vector u are identically zero, since the def- 
inition of u involves only H and not V. Therefore (for cither of these two reasons), 
the amplitudes ak{q) corresponding to the reflection-odd subspace are all identically 
vanishing; these eigenvalues make no contribution whatsoever to the partition func- 
tion. It follows that we can work entirely within the reflection-invariant subspace, 
which has dimension SqCyl(m) [16, Table 2]. For strip widths m > 8, SqCyl(m) is 
strictly smaller than TriCyl(m), so that the matrix T_ is nontrivial. 



4.1 = 2p 

This case is trivial, as the transfer matrix is one-dimensional: 

Z(2p X np) = qiq - l)[(g - 2)(g - 3)]^-' . (4.4) 

Please note that the triangular-lattice strip 2p x np is not equivalent to the strip 
2f X Up for any length n > 2. 

Since there is only one eigenvalue, there is obviously no crossing, hence B = 0. 
However, there are zeros for all n at = 0, 1 and for n > 2 at g = 2, 3. 



4.2 Lx = 3p 

This case is also trivial, as the transfer matrix is again one- dimensional: 

Z(3p X np) = q(q - l)(q - 2)(g^ - V + 29g - 32)^*-^ . (4.5) 

Again B = 0, and the amplitude vanishes at g = 0, 1,2 (which are the first three 
Beraha numbers). For n > 2 there are additional fixed zeros at q ^ 2.5466023485 
and q ^ 3.2266988258 ± 1.4677115087i, where the eigenvalue vanishes. This strip 
was studied by Rocek et al. [41]. 



4.3 Lx = 4p 

The transfer matrix is two-dimensional. In the basis P = {1, S13 + ^24} it can be 
written as 

rrfA ^ - - 12g' + 58^2 - 135? + 126 2(5^ - lOg^ + 33g - 36) \ . . 
^^^^^ - I _2(g2_ 7^ + 12) 2(g2_6g + 9) J ^^"^^ 
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and the partition function is equal to 

Z(4p X np) = q{q - 1) ( ^ ^[/-^ 1^) ^ )^ " ^(^p)""' " ( J ) (4-7) 

The hmiting curve B (see Figure 9) contains three pieces: two complex-conjugate 
arcs and a self-conjugate loop. In addition, at the point q = 3 both eigenvalues 
vanish simultaneously (i.e. the transfer matrix T(4p) itself vanishes); this is a special 
(degenerate) species of isolated hmiting point [27]. (One of the amplitudes does not 
vanish at g = 3, but that is irrelevant.) The self-conjugate loop-like component 
crosses the real axis at go ~ 3.4814056002 and q — 4. There are four endpoints: 

q ^ 1.3705340683 ± 2.7508526144 i (4.8a) 
q ^ 3.6294659317 ± 0.6958422780 i (4.8b) 

This limiting curve was first obtained in the pioneering paper of Beraha and Kahane 
[21] (see also [41]). They drew the important conclusion that g = 4 is a limiting point 
of (complex) chromatic roots for the sequence 4p x np of planar graphs — hence the 
wonderful title of their paper, "Is the Four-Color Conjecture Almost False?" ^ 
The determinant det D{q) has the form 

det D{q) = 8q\q - l)\q - 2)(g' - 3g + l)(g - 3f{q - 4)^ (4.9) 

We recognize the first five minimal polynomials Pk{q) given in [16, Tabic 1]. Hence, 
the determinant vanishes at the first five Beraha numbers g = 0, 1, 2, i?5, 3. It also 
vanishes at g = 4, which corresponds to B^o- The dominant amplitude vanishes only 
at g = 0, 1, 2, B^, these values of g correspond to isolated limiting points. 

The fixed zeros at g = 0, 1 are trivial ones. At the fixed zero q — 2, there is one 
nonzero eigenvalue (A* = 10) with a vanishing amplitude and one zero eigenvalue with 
a nonvanishing amplitude (a = 2); we are therefore in Case 3 described in Section 2. 
The fourth real zero converges exponentially fast to B^ (see Table 3), in agreement 
with the fact that this is an isolated limiting point. The fifth real zero, at g = 3, is a 
fixed zero where both eigenvalues vanish (one of the amplitudes is and the other is 
18); we are therefore on Case 2 described in Section 2. The fact that g = 3 is a fixed 
zero is due to the width not being a multiple of 3. Finally, g = 4 is a crossing point 
where both eigenvalues take the value A = 2 and one of the amplitudes vanishes. 
The sixth real zero in Table 3 converges at an approximate 1/n rate to the value 
go ^ 3.4814056002. 

For this strip there is a vanishing subdominant amplitude at g = i?6 = 3, which is 
greater than B^^i = B^, in contrast with the behavior observed for the square lattice 
[16, 17]. 

^Ironically, by the time that this article was published, the Four-Color Conjecture had become 
the Four-Color Theorem. The Beraha-Kahane article was submitted in 1976 but not published until 
1979. 
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4.4 Lx = 5p 

The transfer matrix is two-dimensional. In the basis P = {1, (^13 + perm.} it can 
be written as 

T(^)-f 5(?^-lV + 76g2-187g+174) \ 

^ ^^P^ - 1^ + iig2 _ 43^ + 58 3^3 _ 35^2 ^ 132^ _ J 

where 

Tn = - 15g^ + 95g^ - 320g^ + 579g - 452 , (4.11) 
and the partition function is equal to 

Z(5pxnp) = 9(?-l)(5-2)(^5[/_%+^^y-T(5pri-( J) (4.12) 

The limiting curve B (see Figure 10) contains three disconnected pieces: two 
complex-conjugate arcs and a self-conjugate loop-like arc. This latter piece crosses the 
real axis at go ~ 3.2072219810 and at g 3.2847747616. There are four endpoints: 

g ^ 0.4772525688 ± 2.5694937945 i (4.13a) 
g 3.5227474312 ± 0.1876729035 i (4.13b) 

The limiting curve was first obtained in [44]. 
The determinant detD(g) has the form 

detL'(g) = 5g'(g - l)'(g - 2)\q^ - 3g + l)(g - 3){q^ - llg' + 43g - 58)' (4.14) 

We recognize the first five minimal polynomials pfc(g) given in [16, Table 1]. Hence, the 
determinant vanishes at the first five Beraha numbers g = 0, 1, 2, B^, 3. Furthermore, 
the dominant amplitude vanishes at all these points; hence they are isolated limiting 
points. The last factor of the determinant does not provide additional isolated limiting 
points. 

The fixed zeros at g = 0, 1, 2 arc trivial ones where all amplitudes vanish due to 
the prefactor g(g — l)(g — 2) in (4.12). The fourth real zero converges exponentially 
fast to (see Table 3). The fifth real zero, at g = 3, is a fixed zero where there is one 
nonzero eigenvalue (A* = —2) with a vanishing amplitude and one zero eigenvalue 
with a nonvanishing amplitude {a = 30); we are therefore in Case 3 described in 
Section 2. The fact that g = 3 is a fixed zero is due to the width not being a multiple 
of 3. Finally, the sixth real zero in Table 3 converges at an approximate 1/n rate to 
the value go ^ 3.2072219810. 



^°Two decades earlier, Beraha, Kahane and Weiss [22] computed the transfer matrix and reported 
the crossings of the limiting curve B with the real axis. But they did not show a plot of the limiting 
curve. 



18 



4.5 Lx = 6p 

The transfer matrix is five-dimensional; it can be found in the Mathematica file 
transfers. m. This strip has been previously studied by Chang and Shrock [44]; but 
they did not compute the limiting curve. 

The limiting curve B is connected (see Figure 11). It crosses the real axis at 
qo pa 3.2524186216. There are four endpoints: 



There are T points at ? 3.9766954928 ± 0.9167681670 i. 
The determinant det D{q) has the form 

det D{q) = 1769472q%q-lf{q-2f{q'^-3q+lf{q-3y\q^-5q^ + 6q-l) 



where the polynomial P{q) can be found in the file transfers. m. The first six 
factors in (4.16) are the first six minimal polynomials given in [16, Table 1]; hence 
the determinant vanishes at the first six Beraha numbers q = 0,1, 2, B^, 3, Bt. It also 
vanishes at the Beraha number Biq, whose minimal polynomial is q^ — bq-\- 5, and 
at g = 4. The dominant amplitude vanishes only at the first six Beraha numbers, 
so these values are the only isolated limiting points. The polynomial P{q) does not 
provide additional isolated limiting points. 

The first two zeros g = 0, 1 are trivial ones. At g = 2 we have three nonzero eigen- 
values with vanishing amplitudes, one zero eigenvalue with zero amplitude, and one 
zero eigenvalue with nonzero amphtude; we therefore fall in Case 3 of Section 2. The 
fourth and fifth real zeros converge exponentially fast to g = S5, 3 (see Table 3). The 
next Beraha number B'j ^ 3.2469796037 is very close to the value go ~ 3.2524186216 
where the limiting curve B crosses the real axis. This explains why the convergence 
rate to the sixth real zero in Table 3 is not as fast as expected (empirically the con- 
vergence is roughly ~ n~^"^); but we expect that it will be ultimately exponential 
(for very large n). We also expect a seventh real zero for large enough n; this zero is 
expected to converge (at an approximate 1/n rate) to the value go ~ 3.2524186216. 
We would need to go to very large n to observe this additional zero. 

Finally, for this strip there is a vanishing subdominant amplitude at g = Bio, 
which is greater than B^^i — B-j, in contrast with the behavior observed for the 
square lattice [16, 17]. 

4.6 Lx = 7p 

The transfer matrix is six-dimensional; it can be found in the Mathematica file 
transf erS.m. 

The limiting curve B is connected (see Figure 12). It crosses the real 
go w 3.4790022937 and g « 3.6798199576. It enters the half-plane Re(g) < 0, and 



g 0.0207708231 ± 2.2756742729 i 
g 4.2838551928 ± 0.6111544521 i 



(4.15a) 
(4.15b) 



x(g2-5g-h5)'(g-4)^P(g)2 



(4.16) 
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there are four endpoints 

q ^ -0.2279183274 ± 2.0134503491 i (4.17a) 
3.9930118897 ± 0.6273386181 i (4.17b) 

There are four T points at g ^ 3.6222949363±0. 1398555812 i and q ^ 3.9816630253± 
0.8993269516 i. FinaUy, it is worth noticing that the hmiting curve encloses a small 
region around 3.479 < Re{q) ^ 3.680 and | lm{q)\ < 0.14. 
The determinant det D(q) has the form 

det D{q) = 6884:U72q^{q - ifiq - 2f{q^ ~ 3q + l^iq - 3f\q^ - 5q'^ + 6q - 1) 
x{q^ -Aq + 2){q^ -bq + bfP{qf (4.18) 

where the polynomial P{q) can be found in the file transfers. m. The first seven 
factors in (4.18) are precisely the first minimal polynomials given in [16, Tabic 1]. 
The next one (g^ — 5g + 5) is the tenth minimal polynomial in [16, Table 1]. Hence, the 
determinant detD{q) vanishes at the Beraha numbers q = 0,1,2, B^,3, B-j, Bg, Biq. 
However, the dominant amplitude vanishes only at the first seven Beraha numbers 
q — 0,1, 2, B5, 3, B7, Bs] these values are the only isolated limiting points for this 
strip. The polynomial P{q) does not provide additional isolated limiting points. 

The first three zeros q = 0,1,2 are trivial ones. The fifth real zero, at g = 3, 
is a fixed zero where there are 2 nonzero eigenvalues with zero amplitudes, one zero 
eigenvalue with zero amplitude, and 3 zero eigenvalues with nonzero amplitudes. The 
fourth, sixth and seventh real zeros converge exponentially fast to q — B^, B^, Bs (see 
Table 3). The eighth real zero seems to converge at an approximate 1/n rate to the 
value go ~ 3.4790022937. 

In contrast with the behavior observed for the square lattice, we find a vanishing 
subdominant amplitude at g = Biq, which is greater than B^^i — Bg. 

4.7 = 8p 

The transfer matrix T(8p) is 15-dimensional. As discussed at the beginning of 
this section, the transfer matrix can be brought into block-diagonal form 

T(8p) ^ ( ^-(«^' ^l^^ ) (4.19) 

where the block T+ (rcsp. T_) is 14-dimensional (rcsp. 1-dimensional) and corresponds 
to the subspacc of rcfiection-invariant (resp. refiection-odd) connectivities. Moreover, 
the amplitude corresponding to the refiection-odd eigenvalue is identically vanishing; 
this eigenvalue therefore makes no contribution to the partition function. The reduced 
transfer matrix T+(8p) can be found in the Mathematica file transf erS .m. The 
one-dimensional block is T_(8p) = — g'^ + 6g^ — 8g — 3. 

There are, however, two further curious features for which we have, as yet, no 
explanation. First of all, we find another eigenvalue A = — g'^ + 6g^ — 8g — 3, this 
time inside the reflection-invariant subspace. Secondly (and even more mysteriously). 



20 



this eigenvalue too has an identically vanishing amplitude. (Wc have checked this fact 
numerically.) The pair of eigenvalues A = — g^ + 6g^ — 8g — 3 can be observed explicitly 
by forming the characteristic polynomial of the transfer matrix, which can be factored 
as^i 

det[r(8p) - Al] = (A + ?'-6?' + 8? + 3)'g2(?,A) (4.20) 

where Q2{q,X) is a polynomial in q and A (it is obviously of degree 13 in A). Un- 
fortunately, we have been unable to find a further change of basis to make T+(8p) 
block-diagonal and thereby bring out explicitly the eigenvalue A = — -|- — 8g — 3 
lying inside that subspace. 

In order to compute the limiting curve B we have mainly used the resultant 
method. The existence of a double eigenvalue A = —q^ + 6q^ — 8g — 3 in the full 
transfer matrix T{Sp) makes the resultant identically zero for ^ = 0. However, this 
problem does not arise if we consider the reduced matrix T+(8p). Nevertheless, the 
existence of an identically zero amplitude within the reduced subspace makes the 
computation of the limiting curve B not completely straightforward, as only those 
eigenvalues with non-identically-vanishing amplitudes should be taken into account 
in computing B. A simple solution is devised by noting that the resultant method 
uses only the characteristic polynomial of the transfer matrix. Therefore, we can 
drop the factor (A + — 6q^ + 8g + 3)^ in (4.20) and compute the resultant using 
the polynomial Q2{q,X)- In this way, we obtain a nonzero resultant, into which the 
zero-amplitude eigenvalues do not enter. 

The hmiting curve B is connected (see Figure 13). It crosses the real axis at 
qo 3.5147694243. It enters the half-plane Re(g) < 0, and there are four endpoints 

q ^ -0.3713655472 ± 1.7983425919 z (4.21a) 
q 4.0496984440 ± 0.7359317819 i (4.21b) 

There are cusp-like structures around q ^ 4.04 ± 0.74 z. A closer look shows that 
these structures are in fact T points located at g ^ 4.0428606488 ± 0.7417105390 i, 
from which there emerge very short branches terminating at the endpoints q ~ 
4.0496984440 ± 0.73593178191 

We form a matrix D{q) of dimension 13 rather than 15; in this way we can avoid 

^^In most of the previously studied cases with cyhndrical boundary conditions (namely, triangular- 
lattice strips of widths 4p < to < 7p and square-lattice strips of widths 4p < m < 9p), the 
characteristic polynomial associated to the transfer matrix cannot be factored as in (4.20). In other 
words, none of the eigenvalues A is a polynomial in q. The cases with to < 3p are trivial as the 
transfer matrix is one-dimensional: there is a single eigenvalue, which is indeed a polynomial in q. 

^^We also checked — though this is not relevant to computing B for the boundary conditions 
being considered here — that the zero-amplitude eigenvalue A = —q'^ + 6q^ — — 3 is not dominant 
at any of the zeros of our resultant. If the zero-amplitude eigenvalue A = —q'^ + Gq"^ — 8q — 3 were in 
fact dominant somewhere in the complex g-plane, then by modifying the top and bottom endgraphs 
(as shown in [41]) it might be possible to make that eigenvalue (in either the reflection-even or 
reflection-odd sector or both) contribute to the chromatic polynomial and thereby obtain a different 
limiting curve B for the different choice of endgraphs. 
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the two identically vanishing amplitudes. Its determinant has the form 



det D{q) = const, x q^^q - lY^q - 2f\q^ - 3g + lf\q - 3)^^^ 

x(g^ - + 6g - lf{q^ - 4g + 2)(g'^ - + 9g - 1) 

x(g2 - 5g + hf\q^ - 7q^ + Uq - 7)\q - 4)^^ 

x(12g^ - 196g^ + 1355g^ - 5126g^ + 11337g^ - 14086? + 7755) 

xPi(q)^ (4.22) 

where the polynomial Pi{q) can be found in the file transfers. m. The first nine 
factors in (4.22) are precisely the first minimal polynomials given in [16, Table 1]; 
therefore, det D{q) vanishes at the Beraha numbers q — B2, . . . , -Bio- The next factor 
is the square of the polynomial q^ — Iq^ + lAq — 7, which is Pi4(g) [16, Table 1], so 
that detD(g) also vanishes 3X q — B14. Finally, the determinant also vanishes at 
q = 4 = Boo- Unlike what we have seen for cylindrical strips of smaller width, in this 
case the remaining part of det D{q) is not the square of a polynomial with integer coef- 
ficients; rather, there is the additional degree-6 factor preceding Pi(g)^. The dominant 
amplitude vanishes only at the first seven Beraha numbers q — 0, 1, 2, S5, 3, S7, Sg, 
so these values are the only isolated limiting points for this strip. The degree-6 factor 
and the polynomial Pi{q) do not provide additional isolated limiting points. 

The first two zeros q = 0,1 are trivial ones. At g = 2, 3, all amplitudes vanish 
except for a few corresponding to zero eigenvalues; we are thus in Case 3 of Section 2. 
At q — 2, there are 10 nonzero eigenvalues with zero amplitudes, 3 zero eigenvalues 
with zero amplitudes, and one zero eigenvalue with nonzero amplitude. At g = 
3 the transfer matrix is not diagonalizable: there are 2 nonzero eigenvalues with 
zero amplitudes, one 2x2 nontrivial Jordan block corresponding to A = with no 
contribution at all to the partition function for all n, and 10 zero eigenvalues with 
nonzero amplitudes. The fourth, sixth and seventh real zeros converge exponentially 
fast to g = BqjBtjBq (see Table 3). The eighth real zero seems to converge at an 
approximate 1/n rate to the value go ~ 3.5147694243. 

Finally, we find vanishing subdominant amplitudes at g = Bio,Bi4, which are 
greater than Bm+i = Bg. 

4.8 = 9p 

The transfer matrix T(9p) is 28-dimensional. It can be brought into block- diagonal 
form 

- ( 't' rie.) ) . (-3) 

where the block T+ (rcsp. T_) is 22-dimensional (resp. 6-dimcnsional) and corresponds 
to the subspace of refiection-invariant (resp. reflection-odd) connectivities. Moreover, 

^•^If we were to form a 15-dimensional or 14-dimensional matrix its determinant would be 

identically zero. This tells us (in case we did not know it already) that two of the amplitudes are 
identically zero. 



22 



all the amplitudes corresponding to the reflection-odd subspacc arc identically van- 
ishing; this subspace therefore makes no contribution to the partition function. The 
reduced transfer matrix T_^{9p) can be found in the Mathematica file transf er3 .m. 

Mysteriously, all of the eigenvalues in the reflection-odd subspace have "copies" 
in the reflection-even subspace. This can be seen by computing the characteristic 
polynomial associated to the transfer matrix T{9p), which can be factored as follows: 

det[T(9p)-Al] = Qi{q,XfQ2{q,X) (4.24) 

where Qi{q,X) and Q2{q,\) are polynomials in q and A; here Qi (resp. Q2) is of 
degree 6 (resp. 16) in A. More specifically, Qi (resp. Q1Q2) is the characteristic 
polynomial of T_ (resp. T+); the fact that Qi appears as a factor in the characteristic 
polynomial of T+ is direct proof of the just-mentioned "copying" of eigenvalues. Even 
more mysteriously, our numerical checks suggest that all the eigenvalues coming from 
Qi have identically zero amplitudes — not only in the reflection-odd subspace (where 
this is well understood) but also in the reflection-invariant subspace. We thus flnd 6 
pairs of equal eigenvalues with identically vanishing amplitudes. 

In order to be able to use the resultant method in this case, we proceed as in 
the previous subsection: we drop the polynomial Qi(g, A)^ (which contains the zero- 
amplitude eigenvalues) from the characteristic polynomial associated to T(9p) and 
compute the resultant with the polynomial (52(5, A) (which contains the eigenvalues 
with nonzero amplitudes) . We have computed the points with ^ = with the resultant 
method; the points with other values of 6 have been computed using the direct-search 
method. 

The limiting curve B is connected (sec Figure 14). It crosses the real axis at 
go ~ 3.5270636990. It enters the half-plane Re(g) < 0, and there are four endpoints 

q ^ -0.4576020413 ± 1.6238415411 i (4.25a) 
q ^ 4.2828643197 ± 0.3823491910 i (4.25b) 

There are T points located at g 4.0160853030 ± 0.78701538591 

We form a matrix D{q) of dimension 16 rather than 28, in order to avoid the 12 
identically zero amplitudes. Its determinant has the form 

det D{q) = const x q^%q - l^^q - 2f\q^ - 3g + l)^-^(g - 3)^^^ 

X (g^ - 5g2 + 6g - lf{q^ - 4g + 2)^(g^ - Qq^ + 9g - 1) 

X (g2 - 5g + 5)^(g^ - 7g^ + 14g - 7fP{q) , (4.26) 

where the polynomial P{q) can be found in the file transf erS.m. We find the same 
"Beraha factors" as in the triangular- lattice strip of width Bp (Section 4.7). Thus, 
det D{q) vanishes at the Beraha numbers g = B2, . . . , -Bio, ^i4- The polynomial P(g) 
is not the square of any polynomial with integer coefficients; rather it can be written 
as Pi{q)P2{qY where Pi and P2 are integer-coefficient polynomials. The dominant 
amplitude vanishes only at the first seven Beraha numbers g = 0, 1,2, 3^,3, B7, Bg] 
thus, these values are the only isolated limiting points for this strip. The polynomial 
P{q) does not provide additional isolated limiting points. 
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The first three zeros q = 0, 1,2 are trivial ones (i.e., all amplitudes vanish iden- 
tically). The fourth, fifth, sixth and seventh real zeros converge exponentially fast 
to q = 3^,3, B-j, Bg (see Table 4). The eighth real zero seems to converge at an 
approximate 1/n rate to the value qo ^ 3.5270636990. 

We again find a vanishing subdominant amplitude at g = B14, which is greater 
than Bm+i = Biq. 



4.9 = lOp 

The transfer matrix 7(10?) is 67-dimensional. It can be brought into block- 
diagonal form 

.(10.) - ( -^(-r) ) . (4.7) 

where the block T+ (resp. T_) is 51-dimensional (resp. 16- dimensional) and corre- 
sponds to the subspace of reflection-invariant (resp. reflection-odd) connectivities. All 
the amplitudes corresponding to the reflection-odd subspace are identically vanishing; 
this subspace therefore makes no contribution to the partition function. The reduced 
transfer matrix T+(10p) can be found in the Mathematica file transf er3 .m. 

The characteristic polynomial of the transfer matrix T(lOp) obviously factors as 
det[T(10p) — Al] = Qi{q, A) Q{q, A), where Qi (resp. Q) is the characteristic polyno- 
mial of T_ (resp. Numerically we have found, once again, that all the eigenvalues 
in the reflection-odd subspace have "copies" in the reflection-even subspace. There- 
fore, the polynomial Q{q, A) should have Qi{q, A) as a factor, so that 

det[T(10p) - Al] = Q,(q, A)' Q2{q, A) (4.28) 

where (^2(9, A) is a polynomial of degree 35 in A. Unfortunately, we have been un- 
able to compute the characteristic polynomial Q{q, A) of the reduced transfer matrix 
T'+(10p) and verify the conjectured factorization (4.28). 

Once again, we have found numerically that the "copied" eigenvalues have identi- 
cally vanishing amplitudes. We thus flnd 16 pairs of equal eigenvalues with identically 
vanishing amplitudes. 

We have used the direct-search method in the computation of B. The hmiting 
curve B is connected (see Figure 15). It crosses the real axis at qo ~ 3.6348299654. 
It enters the half-plane Re(g) < 0, and there are four endpoints 

q ^ -0.510807 ± 1.481233 i (4.29a) 
q 4.113231 ± 0.492835 i (4.29b) 

There are T points located at g ^ 4.0632619066 ± 0.8803786140 i 

We have numerically checked that the dominant amplitude vanishes at the flrst 
nine Beraha numbers q = 0,1,2, B^,3, B-^, Bg, Bg, Biq (and at no others); therefore, 
these values are the only isolated limiting points for this strip. We have found no 
evidence of complex isolated limiting points from the zeros of the finite-length strips. 

The first two zeros q = 0,1 are trivial ones. At g = 2, 3 we are in Case 3 described 
in Section 2. At q — 2, there are 38 nonzero eigenvalues with zero amplitudes, 12 zero 
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eigenvalues with zero amplitudes and one zero eigenvalue with a nonzero amplitude. 
At g = 3 the transfer matrix is not diagonalizable: we find 4 nonzero eigenvalues 
with zero amplitudes, one 2x2 nontrivial Jordan block corresponding to A = —3 
which does not contribute to the partition function for any n, 40 zero eigenvalues 
with nonzero amplitudes, and two nontrivial Jordan blocks (of dimensions 3 and 2 
respectively) corresponding to A = 0. The contribution of these later blocks is zero 
except for n = 1. The fourth, sixth, seventh, eight and ninth real zeros converge 
exponentially fast to g = B^, Bj, Bg, Bg, B^q (see Table 4). The tenth real zero seems 
to converge at an approximate 1/n rate to the value Qq a; 3.6348299654. 

Finally, we have also checked that there are vanishing amplitudes (in addition to 
the trivial 32 identically zero amplitudes) for q = Bn, Bu, and Big. In all these 
cases, the vanishing amplitude corresponds to a subdominant eigenvalue; thus, none 
of these points is an isolated limiting point. Please note that the last two values 
(namely, B14 and Big) are greater than B^^i = Bn. 



The transfer matrix T(llp) is 145-dimensional. It can be brought into block- 
diagonal form 



where the block T+ (resp. TL) is 95-dimensional (resp. 50- dimensional) and corre- 
sponds to the subspace of reflection-invariant (resp. reflection-odd) connectivities. All 
the amplitudes corresponding to the reflection-odd subspace are identically vanishing; 
this subspace therefore makes no contribution to the partition function. The reduced 
transfer matrix T+(llp) can be found in the Mathematica file transf er3 .m. 

The characteristic polynomial of the transfer matrix T(llp) obviously factors as 
det[T(llp) — Al] = Qi{q, A) Q{q, A), where Qi (resp. Q) is the characteristic polyno- 
mial of T_ (resp. T+). Numerically we have found, once again, that all the eigenvalues 
in the refiection-odd subspace have "copies" in the refiection-even subspace. There- 
fore, the polynomial Q{q, A) should have Qi(g, A) as a factor, so that 



where Q2{q, \) is a polynomial of degree 45 in A. Unfortunately, we have been un- 
able to compute the characteristic polynomial Q{q, A) of the reduced transfer matrix 
^^(llp) and verify the conjectured factorization (4.31). 

Once again, we have found numerically that the "copied" eigenvalues have identi- 
cally vanishing amplitudes. We thus find 50 pairs of equal eigenvalues with identically 
vanishing amplitudes. 

As in the preceding subsection, we used the direct-search method in the compu- 
tation of B. This curve crosses the real g-axis at qo 3.6441399017 (see Figure 16). 
It enters the half-plane Re{q) < 0, and there are four endpoints 



4.10 La: = lip 




(4.30) 



det[T(llp)-Al] = Qi(g,A)2Q2(g,A) 



(4.31) 



-0.543988 ± 1.363241 i 



4.156093 ± 0.529420 i 



(4.32a) 
(4.32b) 
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There are T points located at g f» 4.0425923021 ± 0.6927608569 i. 

We have numerically checked that the dominant amplitude vanishes at the first 
nine Beraha numbers q = 0,1,2, B^,3, B-j, Bg, Bg, Biq (and at no others), so that 
these values are the only real isolated limiting points for this strip. We inspected the 
zeros of the finite-length strips and found no evidence of complex isolated limiting 
points. 

The first three zeros q = 0,1,2 are trivial ones, as all amplitude vanish. At 
q — 3 the transfer matrix is not diagonalizable. There are four nontrivial Jordan 
blocks corresponding to four nonzero eigenvalues. Furthermore, the contribution 
of these Jordan blocks to the partition function vanishes for all n > 1. We also 
find 8 nonzero eigenvalues with zero amplitudes, 30 zero eigenvalues with nonzero 
amplitudes, and 30 zero eigenvalues with zero amplitudes. The fourth, sixth, seventh, 
eighth and ninth real zeros converge exponentially fast to g = Br,, Bi, Bg, Bg, Biq (see 
Table 4). The tenth real zero seems to converge at an approximate 1/n rate to the 
value qo ^ 3.6441399017. 

Finally, we have also checked that there are vanishing amplitudes (in addition 
to the trivial 100 identically zero amplitudes) for q = Bu, B12, B14, Big. In all these 
cases, the vanishing amplitude is subdominant; thus, none of these points is an isolated 
limiting point. Again, the values B14 and Bis ^-re greater than B^^i — B12. 



The transfer matrix T(12p) is 368-dimensional. It can be brought into block- 
diagonal form 



where the block (resp. T_) is 232-dimensional (resp. 136-dimensional) and corre- 
sponds to the subspacc of reflection-invariant (resp. reflection-odd) connectivities. All 
the amplitudes corresponding to the reflection-odd subspace are identically vanishing; 
this subspace therefore makes no contribution to the partition function. The reduced 
transfer matrix T+(12p) can be found in the Mathematica file transf er3 .m. 

The characteristic polynomial of the transfer matrix T(12p) obviously factors as 
det[T(12p) — Al] = Qi{q, A) Q{q, A), where Qi (resp. Q) is the characteristic polyno- 
mial of T_ (resp. T+). Numerically we have found, once again, that all the eigenvalues 
in the reflection-odd subspace have "copies" in the reflection-even subspace. There- 
fore, the polynomial Q{q, A) should have Qi{q, A) as a factor, so that 



where Q2{q,X) is a polynomial of degree 96 in A. Unfortunately, we have been un- 
able to compute the characteristic polynomial Q{q, A) of the reduced transfer matrix 
r+(12p) and verify the conjectured factorization (4.34). 

Once again, we have found numerically that the "copied" eigenvalues have iden- 
tically vanishing amplitudes. We thus find 136 pairs of equal eigenvalues with iden- 
tically vanishing amplitudes. 



4.11 = 12p 




(4.33) 



det[r(12p) - Al] = Qi{q,XfQ2{q,X) 



(4.34) 
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Due to the large dimension of the transfer matrix, we have been unable to compute 
the limiting curve. However, we have managed using the direct-search method to 
compute the point where B crosses the real g-axis: go ~ 3.6431658979. We have also 
computed the position of the pair of complex-conjugate T points that are obvious in 
Figure 16; the result is q ^ 4.05713658 ± 0.734324791 

We have numerically checked that the dominant amplitude vanishes at the first 
nine Beraha numbers g = 0, 1, 2, B^, 3, Bj, B^, Bg, and Biq. These values are the only 
isolated limiting points for this strip. We have found no evidence of complex isolated 
limiting points by inspecting the zeros of the finite-length strips. 

The first two zeros g = 0, 1 are trivial ones. At g = 2 there are 164 nonzero 
eigenvalues with zero amplitudes, 2 zero eigenvalues with nonzero amplitudes, and 
66 zero eigenvalues with zero amplitudes. The convergence to g = B^,3,By. Bg, Bq, 
Bio is exponentially fast (see Table 4). The tenth real zero seems to converge at an 
approximate 1/n rate to the value go ~ 3.6431658979. 

Finally, we have also checked that there are vanishing amplitudes (in addition to 
the trivial 272 identically zero amplitudes) for g = Bu, B12, B13, Bu, Big, B22- In all 
these cases, the vanishing amplitude is subdominant; thus, none of these points is an 
isolated limiting point. 

In this case we find three vanishing subdominant amplitudes {Bu, Bis, a-nd B22) 
that are greater than the value B^+i — -B13. 

5 Numerical Results for the Triangular-Lattice 
Chromatic Polynomial: 
Zig— Zag Boundary Conditions 

Until now, we have built up the triangular lattice by transferring along a direction 
that is perpendicular to one of the principal directions of the lattice. An alternative 
choice, of course, would be to transfer along a direction that is parallel to a principal 
direction. When periodic boundary conditions are imposed across the strip, these 
two constructions are inequivalent: they yield different finite graphs. We shall refer 
to this alternative construction, with periodic boundary conditions in the transversal 
direction and free boundary conditions in the longitudinal direction, as "zig-zag" 
boundary conditions, and it will be denoted by the subscript Z. Note that the lattice 
width m must be even. 

For zig-zag boundary conditions, the transfer matrix is not given by the formulae 
of Ref. [16]. Rather, as is evident from Figure 1(b), the transfer matrix now takes 
the following form 

T(mz) = HVevenHVodd , (5.1) 

where Veven (resp. Vodd) is the product of the matrices associated to the vertical bonds 
located at even (rcsp. odd) sites of the lattice. With this definition, all the formulae 
applied in the previous sections hold. 

Our original motivation for introducing this new construction was the following: 
It is clear from Sections 3 and 4 that the limiting curves for the strips with free and 
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cylindrical boundary conditions differ qualitatively by the existence of a small addi- 
tional inward-pointing branch for the case of free b.c, which is absent for cylindrical 
b.c. In the limit of infinite width, one might wonder whether this branch extends 
to g = 2, as the triangular-lattice Ising model is known to have a zero-temperature 
critical point [32] . We found it interesting to examine whether we can recover such a 
branch by imposing zig-zag boundary conditions. The answer turns out to be nega- 
tive; but it seems to us that zig-zag b.c. are interesting in their own right, irrespective 
of this initial motivation. 

We have computed the transfer matrix T(mz) and the limiting curves B for 
triangular- lattice strips of even widths 2 < m = Lj. < 10. It is interesting to note 
that the trick discussed in [16, Sections 3.1 and 3.3] for the standard construction of 
a cylindrical triangular-lattice strip is not necessary here. On the other hand, the di- 
mension of the transfer matrix T(mz) is in general different from TriCyl(m), because 
the invariances are different: for zig-zag boundary conditions, the system is invariant 
under translations of even (but not odd) length and under reflections. 

5.1 Lx = 2z 

This case is trivial, as the transfer matrix is one-dimensional: 

Z(2zXnp)=g(g-l)(g-2)2("-i) (5.2) 

Please note that the strip 2z x is equivalent to 2f x np. Since there is only one 
eigenvalue, there is obviously no crossing, hence B — 0. However, there are zeros for 
all n a,t q — 0,1 and for n > 2 at g = 2. 



5.2 Lx = 4z 

The transfer matrix is three-dimensional. In the basis P = {1, 5i3, 824} it can be 
written as 

(2^11 2^12 Ti3 \ 

2q-5 q^-Aq + A 1 (5.3) 

Tai q^-Aq + A 3?^ - 17g + 25 / 

where 

Tn = g^-12g3 + 56g2-121g+101 (5.4a) 

T12 = - 8q^ + 20q - 16 (5.4b) 

Ti3 = q^ - lOq^ + 34g - 40 (5.4c) 

T31 = 2g3 - 17g^ + 50g - 50 (5.4d) 

and the partition function is equal to 

q^-Sq + sY 

Z{Ap xnz) ^ q{q-l)\ q-l \ ■ T(4z)"-^ • I I (5.5) 
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The limiting curve B (see Figure 18) contains two complex-conjugate arcs, which 
do not cross the real g-axis. There are four endpoints: 

q ^ 2.0991442518 ± 2.5589234827 i (5.6a) 
q fa 2.7371672817 ±0. 1723332852 i (5.6b) 

We have found a complex-conjugate pair of double zeros of the resultant for = 
(see [16, Section 4.1.1]) at q ^ 3.7718445063 ± 1.1151425080 i. At these values the 
limiting curve appears at first glance to be singular (see Figure 18). However, a closer 
look reveals that this is not the case; in fact, the limiting curve is perfectly analytic 
around these two points. -"^^ 

The determinant det D{q) has the form 

det D(q) = -q'(q-mq-2nq'-3q + l)(q-3)\2q-5r 

x(g3- 10^2 + 34^-38)2 (5.7) 

We recognize the first five minimal polynomials pfc(g) given in [16, Table 1]. Hence, the 
determinant vanishes at the first five Beraha numbers g = 0, 1, 2, B^, 3. The dominant 
amplitude vanishes at all these points except at q — 3; therefore, q — 0, 1, 2, are 
isolated limiting points. It is interesting to note that the dominant amplitude vanishes 
also at g = 5/2, so that this too is an isolated limiting point. This is the first time we 
have found a real isolated limiting point that is not a Beraha number. (For the square 
lattice, we did not find any non-Beraha real isolated limiting point [16, 17].) There are 
no additional isolated limiting points coming from the last factor q^ — lOg^ -|- 34g — 38, 
as the vanishing amplitudes in question all correspond to subdominant eigenvalues. 

The first two real zeros g = 0, 1 are trivial as all the amplitudes vanish. The third 
real zero g = 2 belongs to Case 3 of Section 2: the two nonzero eigenvalues have 
zero amplitude, and there is an additional zero eigenvalue with a nonzero amplitude. 
Finally, the fourth and fifth real zeros converge exponentially fast to the values 5/2 
and B5, respectively. In summary, we find five isolated limiting points g = 0, 1,2,5/2 
and i?5- This is in agreement with Table 5. 

Please note that there is a vanishing subdominant amplitude at q = Bq = 3. 
This value is greater than B^^+i — B^, in contrast with the observed behavior for the 
square lattice [16, 17]. 

5.3 Lx = 6z 

The transfer matrix is seven-dimensional; it can be found in the Mathematica 
file transf er3 .m. 



^"^In the notation of [16, Section 4.2], the characteristic polynomial of T{4z) can be expanded as 

P(A, q) = c(A - Ao)2 + d{q - qo f + e{q - go) (A - Aq) + . . . 

around the points qq ~ 3.7718445063 ± 1.1151425080 i and the dominant (double) eigenvalue Ao = 
—0.5237532362 + 4.5580089825 i (the linear terms in A — Ao and q — qo vanish). This expansion leads 
to analytic eigenvalues A±(g) around q = qo and to an analytic equimodular locus around q = qo- 
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The limiting curve B is connected (see Figure 19). It crosses the real axis at 
q ~ 3.1752579126. There are four endpoints: 

q ^ 0.3618461880 ± 2.5093731708 i (5.8a) 
q ^ 4.2589504182 ± 0.7015734543 i (5.8b) 

There are T points at q ^ 3.8395346820 ± 1.1149959335 i. 
The determinant det D{q) has the form 

det D{q) = 81q\q - l)^(g - 2y%q^ -3q+ l)\q - 3)^\q^ - 5q^ + 6q- l)P{qf (5.9) 

where the polynomial P{q) can be found in the file transf er3 .m. The first six polyno- 
mial are the first six minimal polynomials given in [16, Tabic 1] ; hence the determinant 
vanishes at the first six Bcraha numbers g = 0, 1, 2, i?5, 3, B-j. The dominant ampli- 
tude vanishes at the first five Bcraha numbers g = 0, 1, 2, i?5, 3 as well as at three of 
the zeros of P(g), namely q ^ 2.7226328355 and q ^ 3.6696077451 ± 0.95068647361 
This is the first triangular-lattice strip where we find complex isolated limiting points. 
In the square-lattice case, complex isolated limiting points were quite common: we 
found such limiting points for L > 6 with both free and cylindrical boundary condi- 
tions [16, 17]. 

The first two real zeros q = 0,1 are trivial ones. The third real zero q = 2 falls in 
Case 3 of Section 2: there are four nonzero eigenvalues with zero amplitudes, two zero 
eigenvalues with nonzero amplitudes, and one zero eigenvalue with zero amplitude. 
The convergence of the fourth, fifth and sixth real zeros to their corresponding limiting 
values (namely, i?5, 2.7226328355, and 3) is exponentially fast, as shown in Table 5. 
Finally, the seventh real zero converges at an approximate 1/n rate to the value 
g fa 3.1752579126. 

In summary, we find six real isolated limiting points at q — 0, 1, 2, B5, 2.7226328355, 3 
and Bi (see Table 5) and two complex isolated limiting points at q ^ 3.6696077451 ± 
0.9506864736 However, since the complex isolated limiting points are very near B, 
it is very difficult to observe the convergence to them as distinct from the convergence 
to B (see Figure 19). 

5.4 L^ = 8z 

The transfer matrix is 24-dimensional; it can be found in the Mathematica file 
transf erS.m. 

The limiting curve B is connected (see Figure 20). It crosses the real axis at 
q fa 3.3941047539. There are four endpoints: 

q ^ -0.2143469947 ± 2.0301412598 i (5.10a) 
q ^ 4.2899063418 ± 0.5046183096 i (5.10b) 

There are T points at q ^ 4.0055796610 ± 0.8830638824 i. 

In this case we were unable to compute the determinant det D{q). However, we 
checked numerically whether any of the amplitudes vanishes at the Beraha numbers 
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Bn, and if this occurs, whether the vanishing amphtude is the leading one or not. 
We have made this check up to n = 50. We have found that at least one amplitude 
vanishes at the Beraha numbers g = 0, 1, 2, i?5, 3, B^, B8,Bq. The dominant amplitude 
vanishes only at the first six (namely, q = 0, 1, 2, S5, 3, Ey), so that these latter 
numbers are isolated limiting points (see Table 5). In Table 5 we also notice an 
additional isolated zero at g ~ 2.8214204955. We have numerically confirmed that 
this point is indeed an isolated limiting point by minimizing the absolute value of the 
dominant amplitude a* in a neighborhood of that point. 

The first two real zeros g = 0, 1 are trivial ones. At q — 2, there are 13 nonzero 
eigenvalues with zero amplitudes, 2 zero eigenvalues with nonzero amplitudes, and 9 
zero eigenvalues with zero amplitudes. Finally, the other real zeros in Table 5 converge 
exponentially fast to their corresponding limiting values (namely, B^, 2.8214204955. 3, 
and B^). We expect, for sufficiently larger lengths n, an additional real zero larger 
than B7 and converging to q 3.3941047539; but we apparently need to go beyond 
n = 80 to see it. In summary, there are seven real isolated limiting points at q — 0, 
1, 2, B5, 2.8214204955, 3, and B7. 

By minimizing the absolute value of the dominant amplitude a*, we have found a 
pair of complex-conjugate isolated limiting points at g ~ 3. 8327415674±0. 73050211595 i 
(See Figure 20). Again, we are not sure that we have found all complex isolated hm- 
iting points for this strip. 

5.5 = lOz 

The transfer matrix is 87-dimensional; it can be found in the Mathematica file 
transfers. m. 

In this case we were unable to compute the limiting curve. However, we were able 
to compute the value of q where that curve crosses the real axis: q ^ 3.5204366907. 

The matrix D{q) is too large for us to compute its determinant. Instead, we 
have checked numerically the eigenvalues and amplitudes of the transfer matrix at 
the Beraha numbers B^ with 2 < n < 50. We have found that at least one amphtude 
vanishes at the Beraha numbers g = 0, 1, 2, B^, 3, B^, Bg, Bg, Biq, Bu. The dominant 
amplitude vanishes only at the first seven (namely, g = 0, 1, 2, B^, 3, Bj, Bs), so that 
these latter numbers are isolated limiting points (see Table 5). In Table 5 we also 
notice two additional isolated zeros at g 2.8737312493 and g ^ 3.3831285312. We 
have numerically checked that in both cases the leading amplitudes vanish, so they 
too are isolated limiting points. 

The first two real zeros g = 0, 1 arc trivial. At g = 2, there are 35 nonzero eigenval- 
ues with zero amplitudes, 31 zero eigenvalues with zero amplitudes, and 21 zero eigen- 
values with nonzero amplitudes. The convergence of the fourth through ninth real ze- 
ros to their corresponding limiting values (namely, S5, 2.8737312493, 3, Sy, 3.3831285312, B^) 
is exponentially fast as shown in Table 5. Finally, the tenth real zero converges at an 
approximate 1/n rate to the value g ^ 3.5204366907. 

In summary, there are nine real isolated limiting points at g = 0, 1, 2, B^, 
2.8737312493, 3, Sy, 3.3831285312, and Eg. We have been unable to say whether 
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or not there are any complex isolated limiting points; but we do not see any candi- 
date in Figure 21. 



6 Thermodynamic Limit 

In this section we will review the Bethe-Ansatz solution found by Baxter [18, 19] 
for the thermodynamic limit of the zero-temperature triangular-lattice Potts antifer- 
romagnet (Section 6.1), and carefully recalculate Baxter's predictions for the limiting 
curve Boo where the chromatic roots are expected to accumulate (Section 6.2). The 
resulting picture will be substantially similar to that set forth by Baxter [19], but with 
a few important qualitative differences. Next we will compare Baxter's predictions 
for the dominant eigenvalues with our finite-lattice data (Section 6.3), and comment 
on the agreements and discrepancies (Section 6.4). 



6.1 Baxter's solution 

In terms of the variables x and 9 defined by 

q = 2-x-x-^ = 2 + 2cos^ (6.1) 

with |x| < 1 and < Re^ < tt, Baxter defined three functions (eigenvalues) gi{q) as 
follows: 



1 /■°"^,^ sinhfcg f 

log 52(g) = J dk I 



sinhA;^ / sinh[A;(7r-2e)/2] 



(2coshA;6l - 1) sinh(7rA;/2) 

cosh[A;(7r - 2^)/2] 
~ (2 cosh ke + 1) cosh(7rA;/2) 



(6.2b) 



00 



, , , sinh /c^ [sinh /c(7r — ^)] ,^ ^ , 



00 



These formulae were obtained in [18], but the corresponding ranges of validity were 
established only in [19]. In particular, Baxter found that the complex g-plane can be 
divided into three domains X>j [i — 1, 2, 3], in each of which the dominant eigenvalue 
is gi- According to Baxter [19], the intersections of these regions with the real axis 
are as follows: 

Pi n R = > 4} U {g < 0} (6.3a) 
r>2 n R = {qo<q<4:} (6.3b) 
P3 n R = {0<q<qo} (6.3c) 
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Baxter therefore determined the parameter qo by solving the equation 



92{qo) = gsiqo) (6.4) 

via Newton's method and numerical integration of (6.2b)/(6.2c), and found Qq ~ 
3.81967. We have refined this computation using the same method, and find 

?o(Baxter) ^ 3.819671731239719 . (6.5) 

This point is labelled F in Figures 24 and 26 (^f ~ 0.427907971348122). We have 
also plotted the eigenvalues g2 and gs over the entire range < ^ < tt and verified 
that there is only one crossing point, namely (6.5). However, as we will argue later, 
the intersections n R are in fact more complicated than what is claimed in (6.3), 
so that (6.5) is not in fact the correct value of go if Baxter's three eigenvalues are in 
fact the dominant ones. (See Sections 6.3 and 6.4 for further discussion.) 

In order to obtain the limiting curve Boo in the complex g-plane, Baxter [19] took 
advantage of the following simpler expressions for the ratios of eigenvalues: 

- 11 iT^ (6.6b) 



where 



u = e^""'/^ (6.7c) 

In these equations we require that |p| < 1 and \y\ < 1 so that the products converge; 
this corresponds to Im^ < 0. 

As Baxter [19] noted, the parameter 9 enters into the products (6.6a,b) only via 
p or y, respectively, and these two variables are invariant under the transformations 

^^^ + 6A; =^ p^p (6.8a) 

^ +3A;' ^ y^y (6.8b) 



n-9 n-9 

for any integers k,k'. Thus, each solution of \g2/gi\ = 1 in the complex p-plane 
corresponds to an infinite family of solutions in the complex ^-plane: these can be 
thought of as a "primary" solution (namely, the one with largest Re^^ contained in 
the physical region < Re^ < tt) and its "images" under the transformation (6.8a) 
with k > 1. As A; — s> cxD these "image" curves converge to 6* = (g = 4). Likewise, 
each solution of Igs/gi] — 1 in the complex y-plane corresponds to an infinite family 
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of solutions in the complex ^-planc: a "primary" solution (the one with smallest Re 
contained in the physical region) and its "images" under the transformation (6.8b) 
with k' > 1. As A;' — > oo these "image" curves converge to 6' = tt (g = 0). It 
is important to note that the transformations (6.8) do not tell anything about the 
dominant or subdominant character of the equimodular curve at the transformed 
value of ^; this property has to be checked by other means. 

Let us emphasize that the symmetries (6.8) play no essential logical role in our 
analysis: one can, in principle, discover all the equimodular curves by direct search, 
without any reference to these symmetries. But it is useful to know, once one has 
discovered one equimodular curve, that in certain circumstances there must exist 
others as well (and to know their exact location). 

Important Remark. Neither the eigenvalues (6.2b,c) nor the eigenvalue ratios 
(6.6a,b) or Qz/qi = (6.6a)/(6.6b) — nor even their absolute values — are invariant 
under the transformation 9 ^ 9 -\- 27rk" that leaves q invariant. Therefore, we must 
require < Re ^ < tt when using these formulae. 

Remarks. 1. The formula q = 2 — x — maps the disc \x\ < 1 one-to-one onto 
the g-plane cut along the interval [0,4]. Therefore (as Baxter observed [19]), when q 
is real and q > A or q < 0, the parameter x is also real, as is the eigenvalue gi defined 
in (6.2a). 

2. The formula g = 2 + 2 cos 6* maps the strip < Rc^^ < vr one-to-one onto the 
g-plane cut along the intervals (— oo,0] and [4, oo). In particular, when q is real and 
< g < 4, the parameter 9 is real (0 < 6* < tt), as are the eigenvalues g2 and 
defined in (6.2b)/(6.2c). 

3. The definitions (6.2a), (6.6a) and (6.6b) cannot be applied directly on the real 
g-axis for < g < 4; rather, one must consider a limit in which Im^ | and hence 

\p\, \y\ 1 1- 

6.2 Recomputation of Baxter's phase diagram 

In this subsection we will carefully recompute the limiting curve B^o under the 
tentative hypothesis that Baxter's three eigenvalues are in fact the dominant ones. 
(In Section 6.3 we will test this hypothesis against our finite-lattice data.) 

6.2.1 Computation of limiting curves 

Our goal is to compute the locus of points where two or more of the eigenvalues 
gi are equimodular, and to determine at each such point whether these equimodular 
eigenvalues are dominant or subdominant. We carry out this procedure as follows: 

1) We use (6.6a) to compute the locus \g2/gi\ = 1 in the complex p-plane (Fig- 
ure 22); we then transform the resulting plot to the complex ^-plane using 
(6.7a). As noted above, each curve in the p-plane corresponds to an infinite 
family of curves in the ^^-plane. Along each of these latter curves, we compute 
Igs/gil and classify the curve (or portions of it) as dominant or subdominant. 



34 



2) We use (6.6b) to compute the locus \g-i/gi\ = 1 in the complex |/-plane (Fig- 
ure 23); we then transform the resulting plot to the complex ^-plane using 
(6.7b). Each curve in the y-plane corresponds to an infinite family of curves in 
the ^-plane. Along each of these latter curves, we compute \g2/ gi\ and classify 
the curve (or portions of it) as dominant or subdominant. 

3) We use (6.6a,b) to compute the locus \gz/ g2\ = 1 directly in the complex 6'-plane 
and to determine dominance or subdominance. 

4) We combine the three families of equimodular curves into a single 6'-plane plot 
(Figure 24). 

5) Finally, we transform the resulting curves to the complex g-plane using (6.1). 
The resulting "phase diagram" is shown in Figure 25; a detailed view near the 
point g = 4 is shown in Figure 26. 

Despite the explicit formulae (6.6), these computations are far from straightfor- 
ward, due to the slow convergence of the products when \p\ or \y\ is near 1 (i.e. when 
q is near the interval < g < 4 of the real axis) and to the need for very high numer- 
ical precision in intermediate stages of the computation. We discuss these technical 
points in the Appendix. 

The equimodular curves \g2/gi\ = 1 in the complex p-plane are shown in Figure 22. 
Each equimodular curve C„ has exactly two endpoints. As we approach the circle \p\ = 
1, more smaller equimodular curves appear. In order to disentangle the larger curves 
from these new smaller curves, we have followed each equimodular curve carefully as 
it approaches the \p\ = 1 limit. In Figure 22 we have shown all equimodular curves 
that intersect the circle \p\ = 0.99. 

The principal feature is a curve Ci running from point G {p = —1, 6 = ir/G and 
images) via the origin {p = 0, 9 = —iO) to point H (p = —i, 9 = 27r/3 and images). 
The next-longest curve (C2) runs from point I {p — i, 9 — 2t:/9 and images) to point 
J (p = e"''^*/^, 9 = Tn/SO and images). The third-longest curve (C3) runs from point 
K (p = e^'^*/^, 9 = 57r/12 and images) to point L (p = e~'^^^^, 9 = 47r/9 and images). 
Some further equimodular curves and their endpoints are shown on Figure 22 and 
enumerated in Table 8. In this table we have first shown the curves C„ for which both 
endpoints are reasonably well determined (see below) . Then we have listed some other 
well-determined endpoints whose counterparts could not be estimated with sufficient 
accuracy; these points are grouped into the category "Others" . 

It is curious that all these endpoints appear to lie at p = e**^ where is a rational 
multiple of tt (with a small denominator). In order to test this conjecture, we have 
performed detailed fits as follows: For each endpoint we first obtained ten nearby 
equimodular points p with \p\ — 0.990, 0.991, . . . , 0.999. Then we performed a least- 
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squares fit of the data^^ using the polynomial Ansatz 



Argp ArgpFit 



!i + J2<^k{l-\P\)' (6.9) 



TT TT 

fc=l 

in order to estimate the phase (p = Argppn at the endpoint. If the equimodular curve 
is smooth close to \p\ = 1, this Ansatz is expected to work well. We have chosen an 
eighth-order polynomial in order to take into account as many data points as possible 
while allowing a little freedom (we have one degree of freedom in the fits). As a 
check, we have repeated this computation with lower-degree polynomial Ansatze and 
dropping the data with the smallest values of \p\ (in order to have at least one degree 
of freedom in the fit). We have used the stability of the estimates for Aigp^n as a 
guideline to decide whether a fit is good or not (see below) . 

We next asked whether the estimated value of ArgpFit/Tr is or is not close to a 
rational number with a small denominator. We have used the following basic criterion: 
the real number x is "close" to the rational number m/n whenever the 

"discrepancy" = \nx — ml (6.10) 

is sufficiently small. In order to make precise what we mean by "sufficiently small" , 
let us first define 

dnix) — minlnx — ml. (6.11) 

Clearly, for every real number x and every integer n, we have dn{x) < 1/2; so, if we 
are to have good evidence that x = m/n, we will need to insist, at the very least, 
that the "discrepancy" be <^ 1/2. But this requirement is, in fact, nowhere near 
stringent enough: the trouble is that, in our procedure, the denominator n is not 
fixed in advance; rather, we start from the number x and we ask whether there exists 
a (small) integer n such that d„(x) <S 1. And the occurrence of some integers n such 
that dn{x) <^ 1 is by no means a rare event. Indeed, a theorem of Hurwitz (1891) 
asserts that for every real number x, there exist infinitely many integers n such that 
dn{x) < l/{\/5n)}^ So we will need to demand that dn{x) -C 1/n. This suggests 
using the criterion dn{x) < for some 6 > and some s > 1. 

Suppose, in fact, that x is a randomly chosen real number (from a uniform distri- 
bution) . Then it is not hard to show that 

Prob(there exists an integer n with dn{x) < 5/n^) < — -^^-^ (6.12a) 



2^C(g) 

c(^+i; 



(6.12b) 



^^We have made the computations with data truncated to eight decimal digits, hence the error of 
the input data is 10~^. Note, however, that the deviations from (6.9) are not statistical fluctuations; 

rather, they arc "corrections to scaling", i.e. due to neglected higher-order terms in (6.9). A large 
value of the is thus a signal that wc need to include higher-order terms in our Ansatz. 

^^The constant is best possible, as is shown by a; = (1 -|- V5)/2. For a proof of Hurwitz's 
theorem, see [54, Chapter 1], [55, Chapters 5 and 6] or [56, Chapter I]. 
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where 0(n) is the Eulcr toticnt function (i.e. the number of integers 1 < k < n that 
are relatively prime to n) and C(s) is the Riemann zeta function. Therefore, if we 
insist that 

"discrepancy" < '-^^ , (6-13) 

we are following a procedure whose probability of yielding a "false positive signal" 
is less than e (if x is chosen uniformly at random). Otherwise put, we can define a 
"significance level" 

e = l^q^ ^ "discrepancy" . (6.14) 

For concreteness, we have chosen s — 3/2. Roughly speaking, we consider the equality 
X — m/n to be very strongly supported if e ^ 0.001, strongly supported if e ^ 0.01, 
and a plausible guess if e < 0.05. 

We have been able to fit the data corresponding to 42 endpoints. (Wc have nu- 
merically located additional endpoints, but we have not included here those endpoints 
that correspond to very small equimodular curves close to the circle \p\ = 1.) Of these 
42 endpoints, we have obtained a reasonably good fit for the 28 points displayed in 
Table 8: here a "reasonably good fit" is defined as one for which e < 0.05 and for 
which the values of ArgpFit and e are stable under changes in the degree of the poly- 
nomial Ansatz (6.9). We find that 18 points have e < 10~^ (eight of them even have 
e < 10-^), while 10 points have 10"^ < e < 0.05. 

We have found several smaller equimodular curves that are hard to see on Fig- 
ure 22, and we suspect that they too have rational endpoints. Indeed, we would not 
be surprised to learn that such rational endpoints are dense in the unit circle (though 
we have no idea how to prove this conjecture). 

The equimodular curves \g^/gi\ = 1 in the complex y- plane are shown in Fig- 
ure 23. Again, each equimodular curve Dn has exactly two endpoints, and smaller 
equimodular curves appear as wc approach the \y\ = 1 limit. We have shown in 
Figure 22 all equimodular curves that intersect the circle \y\ = 0.99. 

The principal feature is again a curve Di running from point G {y = e^^'^*/^, 
9 — 7r/6 and images) via the origin {p = 0, 6 ~ tt — iO) to point B. {y = 1, 6 = 27r/3 
and images). The next-longest curve {D2) runs from point K (y — e^'^*/^, 9 — 57r/12 
and images) to point L {y = e^'^*/^, 9 = 47r/9 and images). Finally, a small curve {D4) 
nms from point G' {y = e"^^^^^^, 9 = 7r/12 and images) to point H' {y = e~^^'^^/^^, 
9 = 27r/15 and images); we call these points G' and H' because they correspond to 9 

^^Proof of (6.12): Let a; be a uniform random number in [0,1]; we do all our arithmetic in 
the real numbers modulo 1. For any fixed m and n, we clearly have Prob(|na; — rn| < J/n**) = 
min[25/n*+^, 1]. Now we need to sum over pairs (m, n) with 1 < to < n; in doing this, it suffices 
to consider only pairs (TO,n) that are relatively prime, since if m/n = m'/n' with n' < n, then the 
interval \nx — m\ < is entirely contained in the larger interval \n'x — m'\ < bjri!^ . This proves 
(6.12a). The equality with (6.12b) is a standard identity in analytic number theory: see e.g. [57, 
p. 371]. 

Of course, (6.12) is still an overestimate, because it neglects the overlaps between intervals whose 
denominators n and n! do not divide each other. It would be interesting to know what the true 
behavior of this probability is. 
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values that are images of the points G and H under thep-plane transformation (6.8a) 
[though they are not images in the j/-plane]. 

Once again, aU the endpoints appear to he at y = e^^ where ip is a. rational 
multiple of tt with small denominator. Even more curiously, many of these endpoints 
correspond to 9 values that are also observed as endpoints in the \g2/g\\ — 1 plot 
(e.g., all points with a label in Tables 8 and 9). We performed fits to those endpoints 
in the same manner as just explained for the p-plane. We found 44 endpoints, and of 
these we obtained reasonably good fits for the 23 endpoints displayed in Table 9. All 
these points satisfy the criterion e < 0.05, and both Argppit and e are rather stable 
as we vary the degree of the polynomial Ansatz (6.9). Of these 23 endpoints, we find 
12 satisfying e < lO^'' (six of them satisfy the stronger condition e < 10~^), while the 
other 11 endpoints satisfy 10^'^ < e < 0.05. 

We have also found more smaller equimodular curves that are hard to see on 
Figure 23; we conjecture that they have rational endpoints and we wonder whether 
these rational endpoints are dense in the unit circle. 

Remark. It is important to note that the two main points in Figures 24-25 (namely, 
G and H) are extremely well determined in both the p- and |/-planes. In particular, 
their e values range from 4 x 10"^ down to 6 x 10~^. Thus, we can trust that the 
values of 6 for these two points are given exactly by = 7r/6 and 9h — 27r/3 (see 
Tables 8-9). 



6.2.2 Summary of 6'-plane phase diagram 

Let us now describe the resulting zero-temperature "phase diagram" in the com- 
plex ^-plane (Figure 24) and discuss the agreements and discrepancies with respect 
to Baxter [19]. For simplicity we have labelled the points by the same letters as 
in Figure 5 of [19]. The portions of curves where the equimodular eigenvalues are 
dominant (resp. subdominant) are depicted in black (resp. pink). 

The curve A-B (resp. B-C) corresponds to the dominant equimodularity of 
(resp. (72) and gi. These two curves together with the real 9 axis enclose the rest of 
the dominant curves. The position of these points is 

= TT (6.15a) 
9b ^ 0.508588719845180- 0.625516375803391 i (6.15b) 
= (6.15c) 

The point B is triply equimodular (i.e. a T point), so that three equimodular 
curves cross there: 

1) C-B-H, which corresponds to \gi\ = \g2\ (it is dominant along C-B and sub- 
dominant along B-H); and 

2) A-B-Q-G, which corresponds to \gi\ = {g^l (it is dominant along A-B, sub- 
dominant along B-Q, and dominant again along Q-G); 
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3) R-B-Q-F, which corresponds to \g2\ = {gsl (it is subdominant along R-B, 
dominant along B-Q, and subdominant again along Q-F). Point R corresponds 

to = —ioo. 

This last result contradicts [19], where the entire curve B-Q-F is claimed to be 
dominant. The position of these points is^^ 

^ 0.440568708859061- 0.235993788540783 i (6.16a) 
~ 0.427907971348122 (6.16b) 

The point Q is also triply equimodular, so that three equimodular curves cross 
there. Two of them have just been discussed: A-B-Q-G and R-B-Q-F. The third is 
C-Q-G, which corresponds to \gi\ = \g2\', it is dominant along G-Q and subdominant 
along Q-G. Please note that the subdominant curve Q-G lies always to the right of 
(but very close to) the dominant curve Q-G. 

There are infinitely many equimodular curves \gi\ = \g2\ terminating at point G 
{q — 4) and converging to it: they are images under (6.8a) of the two curves G-B-H 
and G-Q-G. The dominant eigenvalue alternates between gi and g2 as these curves 
are crossed. For simplicity, we have shown in Figure 24 only the first few of these 
image curves. The endpoints of these curves can be obtained easily using the values 
of 9g — 7r/6 and 9h — 2tt/3 (see Table 8) and transformation (6.8a). The result is 

= 3(Tfs) 

According to (1.3)/(6.1), they correspond to Beraha numbers: 9G,k is -Bi2+i2fc, and 

0H,k is -B3+i2jfc- 

Likewise, there are infinitely many equimodular curves \gi\ = {g^l terminating at 
point A (g = 0) and converging to it: they are images under (6.8b) of the two curves 
A-B-Q-G and A-H. The dominant eigenvalue alternates between gi and 513 as these 
curves are crossed. Once again, we have shown in Figure 24 only the first few of these 
image curves. The endpoints of these curves can be obtained easily using the values 
of 9g — 7r/6 and 9h — 27r/3 (see Table 9) and transformation (6.8b). The result is 

= ii^. (6.18a) 

In this case, none of them corresponds to a Beraha number. 

There are also many small equimodular curves lying near the real ^-axis, some 
of which are dominant; they arise from the curves I-J, K-L, etc. in Figures 22 and 

^^Dcspitc appearances, the point F in Figure 24 is not triply equimodular: the dominant equimod- 
ular curve l^il = I52I meets the real axis at H' {9 = 27r/15 ~ 0.41888), i.e. slightly below point F; 
this is discussed below at (6.17b) ff. This splitting is somewhat more visible on Figure 26. 
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23. The number of these curves grows rapidly as | Im^| — > 0, so we cannot possibly 
compute all of them; moreover, the computation becomes increasingly difficult as 
|Im6'| — >• 0, since |p| and \y\ are tending to 1. We have, in any case, shown in 
Figure 24 all those curves in the range 0.2 < Re^ < 2.3 that intersect the half- 
plane Im^ < —0.01. Again, one can obtain the endpoints of these curves by applying 
transformations (6.8) to the corresponding 9 values in Tables 8 and 9. These endpoints 
are not in general Beraha numbers. Only point I {9i — 27r/9) and its transformed 
values under (6.8a), 

(6.19) 

' 9 + 12A; ' ^ ^ 

are Beraha numbers (namely, i?9+i2fc)- 

It follows from Figure 24 that the point go should not be identified with F (as 
Baxter [19] did), but rather with the point G at position 

= ^ 0.523598775598299 (6.20) 

(see Tables 8 and 9). 



6.2.3 Summary of g-plane phase diagram 

In Figure 25 we show the above "phase diagram" in the g-plane (for clarity, only 
the dominant cquimodular curves have been depicted). This is quite similar to Fig- 
ure 5 of [19], except for four issues: 

1) The "phase diagram" around point C {q — A) is richer than the one found by 
Baxter (see Figure 26 for a detailed plot of this region). The largest components of 
the region 1^2 (where g2 is dominant) are bounded by C-B-Q-C and its complex- 
conjugate counterpart C-E-Q'-C. The points B and Q take the values 

qB ~ 4.099903170634857 + 0.649694690705481 i (6.21a) 
gq K 3.859627688708099 + 0.203154495450945 i (6.21b) 

However, there are additional components of Vi and 1^2 near point C; indeed, as we 
approach point C the dominant eigenvalue alternates between gi and 5^2 ■ Thus, Qi is 
dominant on the region lying between C-Q-G-Q'-C and C-S15-C; g2 is dominant 
on the region between the curve C--B15-C and the next curve we find towards C, etc. 
Baxter found [via the transformation (6.8a)] only half of the curves around point C, 
namely, those curves that arc images of C-B-H and past through the Beraha numbers 
-B3+i2ifc- (In Figure 26 we have shown two of these curves, corresponding to S15 and 
-627-) 

2) If we define go as the point on the real g-axis where the region I?3 ends, then the 
above discussion implies that go is not given by gp ~ 3.819671731239719 as Baxter 
believed [cf. (6.5)], but rather by 

go = go = S12 = 2 + 3.732050807568877 (6.22) 

[cf. (6.20)]. 
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3) The "phase diagram" around point A (g = 0) is also richer than the one found 
by Baxter. The analytic structure is similar to the one already discussed for the point 
C, except for the fact that now the dominant eigenvalue alternates between gi and g^. 
Again, Baxter found [via (6.8b)] only some of the equimodular curves in this region. 
The physical meaning of these new curves is not clear to us. 

4) We also find many new dominant equimodular curves lying very close to the 
real g-axis between points A and C. Some of them (lying between H and G) have 
been depicted in Figure 25. These curves were missed by Baxter and their physical 
meaning is unclear. 

Let us stress that Baxter was aware of some (though not all) of the "extra" 
equimodular curves discussed here. He said [19, p. 5252]: 

We have ignored those alternative curves as the finite lattice calculations give 
no evidence of there being zeros on them, and the full set of complex zeros just 
fits onto BAE, BCE, BFE. 

It is an open question whether ignoring these curves is indeed the correct thing to do; 
we shall address this question more fully in Sections 6.3 and 6.4. 

On Figure 25 we have also superposed the limiting curve for L = lip (sec Fig- 
ure 16).^^ This curve lies quite close to the main parts of the infinite-strip- width 
limiting curve obtained in this section. The point Qq for the L = lip strip (as well 
as for the rest of the strips considered in this paper: see Table 6) hes below both 
Baxter's prediction (6.5) and our somewhat lower prediction (6.22). So our results 
in Sections 3-5 are compatible with both predictions; unfortunately our strips are 
not yet wide enough to distinguish between them. On the other hand, it is pre- 
cisely around points A and C that the finite-strip limiting curve is not defined, so our 
transfer-matrix results do not give any clue as to whether the additional dominant 
equimodiilar curves we have found there can be neglected (as Baxter did) or not. 

Another way of discovering whether the true value of go is given by point F (6.5) or 
by point G (6.22) is to consider the isolated limiting points of wide triangular-lattice 
strips. We expect that all real isolated limiting points are smaller than qq. Thus, if 
we find any isolated limiting point larger than (6.22) [and smaller than (6.5)], then 
we should conclude that our prediction is false. In Tables 1, 3, 4 and 5 we do not 
find any such zero. The largest real zero we have found is ~ 3.6345747 < qc = -B12 ~ 
3.7320508. . .. On the other hand, for free and cylindrical boundary conditions, all 
real isolated hmiting points are expected to be Beraha numbers. Thus, if go were 
given by F, as Baxter predicts, then the largest real isolated limiting point would be 
Bi4 PS! 3.801938. On the other hand, if go is given by G, as we predict, then the largest 
real isolated limiting point would be Bu ^ 3.682507 (or perhaps B12 ~ 3.732051). In 
Figure 26 we have also depicted the position of the Beraha numbers Bu, . . . , Biq to 
make easier the comparison with the two alternative values of go. 

^^L = 12p would be substantially similar, had we been able to compute the limiting curve for it 
(compare Figures 16 and 17). 
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6.3 Comparison of Baxter's eigenvalues with finite-lattice 
data 

Let us now compare Baxter's predictions for the leading eigenvalues gi{q) [i — 
1,2,3] with our finite-lattice data on strips up to width L = 12 (with cylindrical 
boundary conditions). We shall consider slices of fixed Imq, and plot log as a 

function of Reg. Simultaneously, we shall also plot fk{q,L) — log\Xk{q)\ for the 
five largest (in modulus) eigenvalues of the transfer matrix. Our goal is to see whether 
analogues of Baxter's eigenvalues gi,g2, 93 can be observed in our finite-lattice data. 

Let us begin at large Img, where the behavior is simplest. A plot for Imq = 2 
is shown in Figure 30. Baxter's eigenvalues Qi- 92-93 are shown as a solid blue line, 
a dashed black line, and a dashed-dotted red line, respectively. The eigenvalue gi is 
dominant for Reg ^ —0.52, the eigenvalue 93 is dominant for —0.52 < Reg ^ 3.82, 
and the eigenvalue 91 is once again dominant for Reg ^ 3.82. The eigenvalue 92 
is everywhere strongly subdominant. Comparing these curves with our finite- lattice 
eigenvalues, we see that: 

• In the region Reg ^ —0.52, the leading eigenvalue agrees almost perfectly with 
9i, and the second eigenvalue agrees very closely with g^. 

• In the region —0.52 ^ Reg ^ 3.82, the leading eigenvalue agrees quite well 
(though not perfectly) with g^. Moreover, in a small part of this region near 
its two extremes, we can identify a subleading eigenvalue that appears to cor- 
respond to 9i. However, in most of this region we are unable to identify an 
eigenvalue corresponding to gi — conceivably it exists but is hidden under too 
many other subleading eigenvalues. 

• In the region Reg > 3.82, the leading eigenvalue almost perfectly with gi. 
However, we are unable to identify any eigenvalue corresponding to g^, except 
perhaps for Reg just above the crossover point 3.82. 

We have also extrapolated to the infinite- volume limit L — > 00 the values of the 
free energy fi{L) corresponding to the leading eigenvalue. At values of g where the 
theory is critical, conformal field theory predicts that the finite-size corrections to the 
bulk free energy should be of the form 

ML) = A(oo) + + ... (6.23) 

where c is the system's conformal charge, G is a geometrical factor that equals a/3/2 
on the triangular lattice, and the dots stand for higher-order corrections. At values of 
g where the theory is noncritical, the finite- volume efi'ects should decay exponentially 
as L — > 00, so one may expect that the Ansatz (6.23) will again work well (and simply 
yield a very small number for c). We therefore estimated the free energy /i(oo) by 
inserting the data points L = 11, 12 into the Ansatz (6.23). The results are shown as 
solid green dots on Figure 30. 

For Reg ^ —0.52 and Reg ^ 3.82, the extrapolated free energy /i(cxd) coincides 
almost exactly with the value for L = 12 (as well as with Baxter's eigenvalue gi) , in 
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accordance with our expectation that the region outside the curve ABCEA is noncrit- 
ical and hence has exponentially rapid convergence to the infinite- volume limit. In 
the intermediate regime —0.52 < Reg < 3.82, the extrapolated free energy coincides 
almost exactly with Baxter's eigenvalue g^, (except near the boundaries of this regime, 
where the extrapolation behaves less well); indeed, it coincides with substantially 
better than L = 12 points did, suggesting that the theory in this region is indeed 
critical and satisfies (6.23) with c 7^ 0. 

The behavior is similar at Img = 1.5 (Figure 31) and Imq = 1 (Figure 32); the 
main difference is that the subleading eigenvalue in the region Re g ^ —0.5 is more 
distant from g^ than it was for Imq — 2. In particular, for Img = 1 the two leading 
finite- volume eigenvalues do not even cross; this reflects the fact that Img = 1 lies 
below the leftmost endpoint of the equimodular curve (see Figure 17). However, 
this is almost certainly an effect that will go away for larger strip widths L, as the 
equimodular curve closes up towards q — (compare Figures 28 and 25). The leading 
finite- volume eigenvalue agrees almost perfectly with gi for Req ^ —0.5 and agrees 
quite well with g^ for —0.5 ^ Reg < 4, despite the lack of crossing; moreover, the 
agreement with g^ in the latter interval is almost perfect (except near the endpoints) 
after extrapolation to L = cxd. All this is exactly the behavior one would expect when 
a curve such as Figure 17 is trying to approximate that of Figure 25. The agreement 
in the region Re g ^ —0.5 of the two leading eigenvalues with gi and g^ must therefore 
be reckoned as excellent. Please note that all these slices lie well above the point B 
and the curve AH. 

Consider now the slice Img = 0.5, which crosses the curves BQ and BC and 
twice crosses the curve AH; the corresponding eigenvalues are plotted in Figure 33. 
Here the eigenvalue gi is dominant for Reg ^ —0.56; the eigenvalue g^ becomes 
evcr-so-slightly dominant for —0.56 ^ Reg ^ —0.09; the eigenvalue gi again becomes 
dominant in the interval —0.09 ^ Reg < 0.95 (between the two crossings of the curve 
AH), where it has a sharp peak; the eigenvalue g^ is dominant for 0.95 ^ Reg ^ 4.01; 
the eigenvalue g2 is briefly dominant in the interval 4.01 < Reg ^ 4.17; and flnally, 
the eigenvalue gi is dominant for Reg ^ 4.17. But very little of this complexity is 
seen in our flnite-lattice eigenvalues: 

• In the region Reg < —0.56, the leading eigenvalue agrees almost perfectly with 

• In the entire region —0.56 < Reg < 4.01, we see no sign whatsoever of gi; in 
particular, we see no sign of the sharp peak in the interval —0.09 ^ Reg ^ 0.95. 
Rather, the leading eigenvalue everywhere agrees closely (though not perfectly) 
with gs. 

• The region Re g fa 4 is very hard to disentangle. 

• For Reg > 4.17, the leading eigenvalue agrees almost perfectly with gi. There 
may conceivably be one or more eigenvalues corresponding to g2', this is uncer- 
tain. 
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When we extrapolate the leading eigenvalne to L — > cxo, there is essentially no 
change in the regions Reg < —0.56 and Reg > 4.17 (where the agreement with 
gi was essentially perfect anyway); on the other hand, we obtain an even better 
agreement with in the intermediate region —0.56 < Reg < 4.01, except near the 
endpoint Re g ~ 4 where the extrapolation goes haywire. 

The behavior is similar at Img = 0.35 (Figure 34). The peak in gi is even sharper, 
and once again no sign of it is seen in the finite-lattice data. Please note that the 
slice at Img = 0.35 intersects only the outermost curve AH, but it is very close to the 
next curve AH' (the maximum value of Img along this second curve is ~ 0.3435). 

Likewise at Img = 0.2 (Figure 35): this slice intersects two of the curves (AH 
and AH'), corresponding to a sharp dip in gi near Reg ~ 0.36, followed by an even 
sharper rise near Reg f« 0.89 and then another precipitous dip. Once again, we see 
no sign of this behavior in the finite-lattice eigenvalues or in the extrapolated leading 
eigenvalue. The same situation holds for Img = 0.1 (Figure 36): there are more and 
sharper dips and rises in gi, including some near Reg ~ 4; but there is no sign of this 
abrupt behavior in the finite-size data or the extrapolated data. 

Finally, let us consider Img = (Figure 37). Here gi is well-defined only for g < 
and g > 4 [cf. (6.2a)], while g2 and g^ are well-defined only for < g < 4 [cf. (6.2b,c)]. 
The eigenvalues g2 and gz cross at the point gp ~ 3.819671731239719 found by Baxter 
[cf. (6.5)]. The leading eigenvalue /i for L — 12p reproduces well Baxter's prediction 
for the leading eigenvalue nearly everywhere (namely, gi for g < and q > 4:, g^ for 
< g < gr, and g2 for Qf < q < 4). The main discrepancy occurs around g = q-p: we 
find that our finite-size data agree with g2 in the interval qo{L = 12p) pa 3.64317 < 
g < 4, which is somewhat larger than the interval gp ^ 3.81967 < g < 4 predicted by 
Baxter. But this is perfectly understandable as a finite-size effect due to the fact that 
qo{L = 12p) is smaller than gp. If we look at the extrapolation to infinite volume of 
the leading eigenvalue, we observe an even better agreement with Baxter's prediction 
in the interval ^ Reg ^ 3; but in a wide interval around g = gp the extrapolation 
looks horrible. This is probably due to the large finite-size effects that characterize 
that region where three eigenvalues are closely competing for dominance and there is 
a complicated pattern of level crossings. 

In summary, our finite-lattice data are in excellent agreement with Baxter's predic- 
tions, except on one major point: we find no evidence of an eigenvalue corresponding 
to gi in the region of the g-plane enclosed by the curve ABCEA. As a consequence, 
we find no evidence of the equimodular curves in this region that should correspond 
to the crossings of gi and g^ (e.g. AHA and its images, LK, etc.) or of gi and g2 (e.g. 
CQGQ'C and its images, etc.). There seem to be two possibilities: 

(a) The eigenvalue gi really is present in this region, but only for strip widths L 
larger than those we have studied (possibly much larger). In this case, the 
limiting curve really would exhibit all the complexities shown in Figures 25 
and 26, and the correct value of go would be given by gc = 2 -|- ^/3 3.73205 
rather than by gp ^ 3.81967 as predicted by Baxter. 

(b) For some reason, the eigenvalue gi is not present in this region (though it is 
clearly present elsewhere) . In this case, the limiting curve Boo would be given by 
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ABCEA and CBQFQ'EC only, and the correct value of go would be qf ^ 3.81967 
after all. 

At present our numerical data are insufficient to distinguish definitively between these 
two scenarios. But some relevant theoretical considerations will be presented in the 
next subsection. 

6.4 Critical discussion 

Baxter [19, p. 5242] stresses that 

we expect the results to be exact: the only way they could be wrong would be if 

the domain structure were more compUcated. For instance, we cannot rule out 
the existence of a fourth domain within which W{q) [the "partition function 
per site"] has some yet different form. Ah we can say is that we have seen no 
sign of it. 

But the trouble seems to be not that a fourth dominant eigenvalue appears in some 
part of the complex plane — we see no sign of that either — but rather that in part 
of the complex plane one of Baxter's three eigenvalues (namely, gi) does not appear! 

One possible explanation may have been hinted at by Baxter in his first paper, 
where he points out [18, p. 2823] that the equivalence between Nienhuis' loop model 
and the triangular- lattice chromatic polynomial "should be treated with caution"; 
indeed he observes that it gives manifestly wrong answers on a finite lattice when 
q — 1 or 3. He then gives the following explanation for this discrepancy, which we 
now find strikingly prescient: 

We can resolve this apparent contradiction by repeating Nienhuis' argument 
for a finite lattice. We start by writing the triangular Potts model as a Kagome 
lattice six- vertex model, as has been done by Baxter et al (1976).^*^ Wc inter- 
pret this as an SOS model by placing spins on the faces of the Kagome lattice, 
adjacent spins differing by 5, the greater being to the left of the intervening 
arrow. Choosing the hexagonal (triangular) faces to have integer (half-intcgcr) 
spins and taking J = —00, we can sum over the half- integer spins to obtain 
Nienhuis' SOS form of the loop model, hut with special boundary weights [em- 
phasis added). For Q < q < A these boundary weights are complex, so it is 
possible for them to modify the bulk behaviour. Obviously this is happening 
in the q = 1 and 3 cases mentioned above. 

In terms of the loop model transfer matrices, these boundary conditions 
may mean that not all eigenvalues contribute to the partition function. (If 
an eigenvector is 'orthogonal to the boundary vector', then it never enters the 
calculation of Zloop-) This must be happening in the q = 1 and 3 cases, so is 
presumably a general phenomenon. (Similar problems with using the six- vertex 
form of the Potts model have previous[ly] been observed (see Baxter 1982b). 

^"The reference being cited hero is Baxter, Kelland and Wu [58]. 
^^The reference being cited here is [59]. 
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Thus not only is care necessary to select the largest eigenvalue of the loop 
model transfer matrix but one should also verify that it is a contributing eigen- 
value. 

Indeed, it is quite possible that Baxter's solution of the hexagonal-lattice loop 
model corresponds to some unusual boundary condition for the triangular-lattice 
chromatic polynomial — different from the free and cylindrical boundary conditions 
studied here — in which one or more eigenvalues present in the latter are entirely 
absent in the former (or vice versa). The presence or absence of such eigenvalues 
can, of course, radically change the limiting curve B in case the "missing" eigenvalue 
becomes dominant somewhere in the complex g-plane.^^ 

But this explanation has a severe defect: if true, it would suggest that the "miss- 
ing" eigenvalue should be absent throughout the complex q'-plane (corresponding to 
an identically vanishing amplitude). In our case, by contrast, the eigenvalue gi is 
unambiguously observed outside the curve ABCEA; it is only inside this curve that 
the eigenvalue gi "disappears" . It is hard to see how such a g-dependent effect could 
be caused by boundary conditions. 

There remains, of course, the other possibility noted earlier: namely, that the 
eigenvalue gi really is present inside the curve ABCEA as well, but only for strip 
widths L larger than those we have studied. Indeed, it is not so surprising that 
on lattices L < 12 we have failed to see the gyrations of gi that give rise to the 
equimodular curves AH, AH', . . . : after all, in this same neighborhood {q ~ 0) we 
also fail to see the crossing of gi and gs that corresponds to the completion of the 
equimodular curve BA at q — 0, and there is little doubt that this feature will be 
restored in the limit L oo. 

Light could perhaps be shed on this issue by a more detailed examination of 
Baxter's Bethe Ansatz for finite- width strips [18] and a comparison of his finite- width 
results with our own. It could also be useful to study boundary conditions that are 
periodic in the longitudinal direction (e.g. cyclic or toroidal), as these seem to lead 
to curves B that close at g = already for finite L [45, 46, 47]. The trouble is that 
these studies are restricted, at present, to very small L (namely, L < 4). 

^^Here is an example of such a situation: In computing the triangular-lattice chromatic polynomial 
with cylindrical boundary conditions, one has the choice of working entirely within the translation- 
invariant subspace (as we have done) or else working in a larger space that includes non-translation- 
invariant connectivities (as Rocc;k et al. [41] have done). If one takes the latter approach, one will 
find that one of the eigenvalues has an identically vanishing amplitude (provided that translation- 
invariant endgraphs are used); when this non-contributing eigenvalue is ignored (as it should be), 
both approaches will give the same limiting curve B. On the other hand, the latter approach 
allows one to consider also the use of non-translation- invariant endgraphs; and it can sometimes 
be arranged for the "extra" eigenvalue to contribute and indeed be dominant in some region of the 
complex g-plane, with the result that the limiting curve B is different from what it was without that 
eigenvalue. Compare, for example, [41, Figures 4(a) and (b)], which show the limiting curves B for 
the triangular lattice of width L = 4p with different endgraphs. 
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7 Summary and Outlook 



7.1 Behavior of dominant-eigenvalue-crossing curves B 

In this paper we have computed the transfer matrix for triangular-lattice strips 
of width 3 < < 9 for free boundary conditions, 4 < L^- < 12 for cylindrical 
boundary conditions, and = 4, 6, 8, 10 for "zig-zag" boundary conditions. The 
transfer matrix allows the computation of the chromatic zeros for strips of arbitrary 
length Ly. As the length Ly tends to infinity (for fixed width Lj.), the chromatic 
zeros accumulate along certain curves (limiting curves B) and around certain points 
(isolated limiting points) according to the Beraha-Kahane- Weiss theorem [25, 26, 27]. 
For all the above strips except = IOf, 12p, lOz, we have been able to compute the 
limiting curves B. The exact computation of all the isolated limiting points has been 
carried out for < 6f, < 9p and < 6z', for the larger strips we were able to 
check that certain values of q arc isolated limiting points, but we cannot be certain 
that we have found all of them. By studying the behavior of the limiting curves and 
isolated limiting points as a function of the strip width (and boundary conditions), 
we hope to shed light on the thermodynamic limit L^iLy — > oo. 

The basic properties of both limiting curves and isolated limiting points are sum- 
marized in Table 6. In all cases the identity 

endpoints = (2 x components) + (2 x double points) + (T points) 



holds. This identity can be derived by simple topological/graph-theoretic arguments. 

By inspection of Table 6, we observe some regularities when becomes large. 
For all three boundary conditions, the curve B appears to become connected (#C = 
1) when is large enough. For all L^, the number of endpoints (#E) is 6 for free 
boundary conditions and 4 for the other two boundary conditions. (Note, however, 
that for Lx > 7f, > lOp and = 8z our counts on the number of endpoints are 
only lower bounds: we may have missed some.) When is large enough, the number 
of T points (#T) is 4 for free boundary conditions and 2 for the other boundary con- 
ditions. We have found no evidence of double points for any of the strips considered. 
Finally, the number of enclosed regions is zero except for — 5^, < 5p and 
Lx = 7p. These regularities are in sharp contrast with the square-lattice case [16, 17], 
where the number of connected components and endpoints seems to grow with the 
strip width. In particular, we have not found in the triangular-lattice limiting curves 
any trace of the small gaps and bulb-like regions that are so common in the square- 
lattice case. It therefore seems that the thermodynamic limit may be achieved in a 
smoother way for the triangular lattice than for the square lattice. Finally, it is worth 
mentioning that in all cases except = 4f and 4z, the limiting curve B crosses the 
real g-axis, thus defining qq. By contrast, for the square lattice, Qq is well-defined only 
for odd widths; for even widths with both free and cylindrical boundary conditions, 
we found either that B fails to intersect the real axis or that it contains a segment of 
the real axis passing through a double point [16, 17]. 




(7.1) 
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Let us also note that, as in the square-lattice case [16, 17], we find chromatic zeros 
with Reg < 0. Indeed, for L^. > 7p and L^. > 8z we find that the limiting curve B 
intersects the half-plane Reg < 0. For free boundary conditions, none of our limiting 
curves {L^ < Sp) reach this half-plane; but from Table 6 we can see that minReg is 
decreasing and we expect that it will be < for > lOp (and possibly already for 
9f). 

The regularities exhibited by the limiting curves become clearer when we super- 
pose them all (with fixed boundary conditions). This is done in Figures 27, 28 and 29 
for free, cylindrical and zig-zag b.c, respectively. We find an overall behavior similar 
to that found for the square lattice [16, 17]. For free boundary conditions (Figure 27), 
we find a monotonic behavior with the width L^: both the leftmost arcs and the right- 
most arcs move outwards as is increased (see also the columns labelled minReg 
and max Re g in Table 6). The value of go (or Re go for = 4f) is also monotonically 
increasing in L^. The overall shape of the limiting curves is similar to the expected 
limiting curve in the thermodynamic limit (Figure 25). We expect that as grows, 
the leftmost endpoints will tend towards g = 0, while the rightmost endpoints will 
go to g = gc(tri) = 4. The crossing point go will eventually go to either point F [cf. 
(6.5)] or point G [cf. (6.22)] in Figure 25. Unfortunately, our numerical data are not 
good enough to tell unambiguously the true limit. There is additional one feature of 
the limiting curves with free boundary conditions that does not correspond to any 
feature of the predicted thermodynamic-limit curve: namely, a pair of small complex- 
conjugate branches emerging from T points and pointing inwards. From Figure 27, 
it seems that the size of these branches does not go to zero as is increased (at 
least up to 8 or 9); rather their size stays more or less constant. We are unable to 
say whether these branches will get shorter for larger and ultimately disappear in 
the limit L^^ — > oo. 

In Figure 28 we superpose all the limiting curves with cylindrical boundary condi- 
tions. As in the square-lattice case, the behavior of the leftmost part of these curves 
seems to be monotonic: the arcs move outwards increased. In particular, 

minReg decreases monotonically with the strip width (see Table 6). However, the 
behavior on the right side of the plot is clearly not monotonic: there are differences 
depending on the quantity mod 3. This is to be expected, since with periodic 
boundary conditions in the transversal direction, strip widths that are not multiples 
of 3 are somewhat unnatural as they introduce frustration in the 3-state Potts an- 
tiferromagnet. Thus, the dependence on in the interval 3 ^ Reg ^ 4.5 is not 
a surprise (the same feature is present in the square-lattice case, where we find an 
even-odd dependence on the limiting curves [16, 17]). For fixed values of mod 3, 
we find that go is monotonic in L.^-: for = 1 mod 3 it decreases, while for = 
or 2 mod 3 it increases. The shape of the limiting curves for = mod 3 is very 
similar to the expected thermodynamic limit (Figure 25), while in the other cases 
we find stronger finite-size effects that we expect to disappear in the limit ^ oo. 
Finally, it is worth mentioning the absence of the extra branches that appear for free 
boundary conditions. 

In Figure 29 we superpose the limiting curves for "zig-zag" boundary conditions. 
In this case we only have three curves, so we are unable to extract any definitive 
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conclusion. We can only confirm the monotonic behavior on the leftmost side of the 
curves and the absence of any extra branches. Again, the overall shape is similar to 
the expected thermodynamic limit depicted in Figure 25. 

If we compare the limiting curves for different boundary conditions, we see that 
the thermodynamic limit is achieved faster for cylindrical and "zig-zag" boundary 
conditions than for free boundary conditions (due to existence of "surface" effects in 
the later). This is similar to the behavior observed in the square-lattice case. 

The main unsolved problem is, of course, whether the extra equimodular curves 
predicted by Baxter's formulae (Figure 25) are really there for large enough L. Un- 
fortunately, the evidence from strip widths L < 12 is inconclusive (Sections 6.3 and 
6.4). 

Remcirk. In the computation of the limiting curves B for triangular-lattice strips 
with cylindrical boundary conditions, we have found a curious behavior involving 

identically vanishing amplitudes. (This is important, because the eigenvalues corre- 
sponding to the identically vanishing amplitudes must be excluded from the compu- 
tation of the equimodular curves. For square-lattice strips, by contrast, we have not 
observed any identically vanishing amplitudes [16, 17].) As explained in the intro- 
duction to Section 4, the transfer matrix can be written (after a change of basis) in 
the block-diagonal form 

T(-) - ( '^r' T.(ip) ) • p-^' 

where the matrix T+(mp) lives on the subspacc of reflection-invariant connectivi- 
ties and has dimension SqCyl(m) [i.e., the dimension of the transfer matrix for a 
square-lattice strip of width m with cylindrical boundary conditions], while the ma- 
trix T_(mp) lives on the subspace of reflection-odd connectivities and has dimension 
TriCyl(m) - SqCyl(m). For m > 8p we have TriCyl(m) > SqCyl(m) and this de- 
composition becomes nontrivial. Now, simple symmetry arguments (see Section 4) 
explain why all the eigenvalues in the reflection-odd subspace should have identically 
vanishing amplitudes. But what is curious and mysterious is that each of these eigen- 
values has an identical "partner" in the reflection-even subspace, also with identically 
vanishing amplitude. This means that the characteristic polynomial associated to the 
transfer matrix T(mp) can be factored as 

det[T(mp)-Al] = Q^iq, Xf Q^iq, X) , (7.3) 

where the zero-amplitude eigenvalues are those coming from the factor Qi{q, A)^. In 
particular, the number of eigenvalues with zero amplitude (# VA) is always even, 
and it equals twice the dimension of the reflection-odd subspace: 

#VA(m) = 2 X [TriCyl(m) - SqCyl(m)] (7.4) 

(see Table 7). This is, at any rate, what we have found for 8p < m < 12p (see 
Sections 4.7-4.11); we conjecture that it holds for larger widths as well. It follows 
that the effective dimension of the transfer matrix TriCyl'(m) is given by 

TriCyl'(m) = 2SqCyl(m) - TriCyl(m) . (7.5) 
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Numerical values for all these quantities are displayed in Table 7; of course, the values 
of 7^ VA and TriCyl'(m) for m = 13, 14 are conjectures. 

Numerical values of TriCyl(m) and SqCyl(m) were first reported in [16, Table 2]. 
An analytic formula of TriCyl(m) for prime values of m has been obtained in [53, 
Theorem 3]. This paper also contains a conjecture for SqCyl(m) with prime m > 3 
[53, Conjecture 2]. Finally, an analytic formula for TriCyl'(m) has been conjectured 
for arbitrary values of m [53, Conjecture 1]. 

7.2 Behavior of amplitudes and the Beraha conjecture 

Let us now discuss the isolated limiting points and the role of the Beraha numbers 
in the triangular-lattice strips. Our results show that the number of isolated limiting 
points is a non- decreasing function of the strip width (for each boundary condi- 
tion), at least up to the maximum we have been able to investigate. For free and 
cylindrical boundary conditions we did not find any complex isolated limiting points 
(see Table 6). For "zig-zag" boundary conditions we find a pair of complex- conjugate 
isolated limiting points for = 6z, and we have evidence of the existence of another 
pair of complex-conjugate isolated limiting points for = Sz- 

Concerning the real isolated limiting points, most of them are Beraha numbers 
(1.3). It is only for "zig-zag" boundary conditions that we find real isolated limiting 
points that are not Beraha numbers: for Lx = 4z we find q = 5/2; for = 6z, 
q ^ 2.722633; for = 8z, q ^ 2.821420; and for = lOz, we find two such points, 
q pa 2.873731 and q ~ 3.383129. It is not clear to us how these non-Beraha real 
isolated limiting points behave as — > oo, e.g., whether their number is bounded or 
unbounded. 

For all the lattices we have studied, we observed empirically that there is at 
least one vanishing amplitude ai{q) at each of the Beraha numbers up to (see 
Table 10). It is reasonable to conjecture that this holds for all L (in agreement with 
a similar conjecture for the square lattice [16, Conjecture 7.1]): 

Conjecture 7.1 For a triangular-lattice strip of width L with free, cylindrical or 
"zig-zag" boundary conditions, at each Beraha number q = B2, ■ ■ ■ , Bl+i there is at 
least one vanishing amplitude ai{q). That is, det D{q) = for q = B2, . . . , Bl+i- 

In contrast with the square-lattice case [16, Conjectures 7.2 and 7.3], however, we 
find that there is a vanishing amplitude [hence det D{q) = 0] also at some Beraha 
numbers larger than Indeed, we find examples for each boundary condition 

(see Table 10): 

• For free boundary conditions, q — Bq is a, zero of det D{q) for L — Ap. 

• For cylindrical boundary conditions, g = is a zero of det D{q) ior L — ip; q — 
Bio is a zero for L — 6p, 7p, Bp; q — Bu is a zero for L — Bp, 9p, lOp, lip, 12p; 
q — B18 is a zero for L — lOp, lip, 12p; and finally, q — B22 is a zero for L — 12p. 

• For "zig-zag" boundary conditions, q — Bq is a, zero of det D{q) for L = 4z. 
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We have systematically checked all Beraha numbers up to B50 to make this list. 

Please note that for free and "zig-zag" boundary conditions, we have found only 
one case each in which a Beraha number beyond -Bl+i is a zero of det D{q), namely 
the relatively small value of L = 4. It is conceivable that for all larger L there are 
no such Beraha zeros, as is conjectured for all L for the square lattice [16, Conjec- 
ture 7.2]. However, this conjecture clearly cannot be true for cylindrical boundary 
conditions. There is presumably some pattern that tells us which Beraha numbers 
q = Bk with k > L + 1 can be zeros of detD[q). Thus far only a few Beraha numbers 
(i?6, -Bio, -B14, -B18, -B22) have appeared on that list. Indeed, we conjecture that the 
pattern is the following: 

Conjecture 7.2 For a triangular-lattice strip of width L with cylindrical boundary 
conditions, the Beraha numbers where detD(g) vanishes is given by the union of the 
sets {^2, S3, S^+i} and {B4^k_2 \ k — 1,2,..., [L/2J}, the upper limit on k being 
the integer part of L/2. 

7.3 Nature of the fixed zeros 

In Section 2 we discussed the "fixed" zeros that occur at small integers q (here 
g = 0, 1, 2, 3) when the graph fails to be g-colorable. From the point of view of the 
transfer-matrix formalism, these fixed zeros can arise in either of three ways: 

1) All the amphtudes vanish at q. Then — for all lengths n>l. 

2) All the eigenvalues Afc vanish at q. Then = for all n > 2. 

3) "Mixed case": Neither all the amplitudes nor all the eigenvalues vanish at q, 
but for each k either or vanishes at q (or both). Then Z^^q) = for all 
n>2. 

Let us now summarize what we have found concerning the nature of these fixed zeros 
for triangular-lattice strips: 

g = 0, 1. At g = 0, 1 all the amplitudes vanish, due to the prefactor q{q — 1) in 
the left vector u. These points therefore belong to Case 1. 

q = 2. At q = 2 the behavior depends on the boundary conditions and on the 
strip width L^: 

• Free boundary conditions: 

— Lj; — 2f: The one eigenvalue vanishes at q — 2 (Case 2). 

— Lx — 3f: There is one nonzero eigenvalue with a zero amplitude, and one 
zero eigenvalue with a nonzero amplitude (Case 3). 

— = 4^: There is at least one nonzero eigenvalue with a zero amplitude 
and exactly one zero eigenvalue with a nonzero amplitude (Case 3). The 
transfer matrix at g = 2 is not diagonalizable: it has a nontrivial Jordan 
block corresponding to A = 0. 
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— Lx > Bp: There is at least one nonzero eigenvalue with a zero amplitude, 
at least one zero eigenvalue with a zero amplitude, and exactly one zero 
eigenvalue with a nonzero amplitude (Case 3). We also find that for all 
Lx > 4f the transfer matrix at g = 2 is not diagonalizable, i.e. it has 
nontrivial Jordan blocks (all corresponding to eigenvalue A = 0). 

• Cylindrical boundary conditions: 

— Lx odd: All the amplitudes vanish, due to the prefactor q{q — l){q — 2) in 
the left vector u (Case 1). 

— Lx — 2p: The one eigenvalue vanishes at q = 2 (Case 2). 

— = 4p: There is one nonzero eigenvalue with a zero amplitude, and one 
zero eigenvalue with a nonzero amplitude (Case 3). 

— Lx even > 6p: There is at least one nonzero eigenvalue with a zero ampli- 
tude, at least one zero eigenvalue with a zero amplitude, and at least one 
zero eigenvalue with a nonzero amplitude (Case 3). 

• Zig-zag boundary conditions: 

— Lj. — 2z: This is identical to — 2f (Case 2). 

— L-r = 4z: There are two nonzero eigenvalues with zero amplitudes, and one 
zero eigenvalue with a nonzero amplitude (Case 3). 

— Lx even > 6z' There is at least one nonzero eigenvalue with a zero ampli- 
tude, at least one zero eigenvalue with a zero amplitude, and at least one 
zero eigenvalue with a nonzero amplitude (Case 3). 

q = 3. The point g = 3 is a fixed zero only for cylindrical boundary conditions 
with strip widths that are not a multiple of 3. There are two distinct situations: 

— = 4p: Both eigenvalues vanish, so that the whole transfer matrix van- 
ishes (Case 2). 

— Lx = 5p: There is one nonzero eigenvalue with a zero amplitude, and one 
zero eigenvalue with a nonzero amplitude (Case 3). 

— Lx = 7p,8p, lOp, lip: There is at least one nonzero eigenvalue with zero 
amplitude, at least one zero eigenvalue with a nonzero amplitude, and at 
least one zero eigenvalue with a zero amplitude [or nontrivial Jordan block 
corresponding to A = with no contribution to the partition function for 
any n > 1] (Case 3). We also find that for > 8p the transfer matrix at 
g = 3 is not diagonalizable, i.e. it has nontrivial Jordan blocks. 
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oo 

A Numerical Computation of 11(1"" tx^) 

n=l 

In this appendix we discuss briefly some of the technical issues involved in the 
numerical computation of Baxter's products (6.6). Everything can be expressed in 
terms of the function 

oo 

R{t,x) = n(l-^^"), (A.l) 

n=l 

which is defined for complex t and x satisfying |a;| < 1 and was first studied by Euler 
[60]. Here we need the cases t — ±1 and t — ± a, cube root of unity. A more detailed 
discussion, including proofs, can be found in [61]. 

The numerical computation of R{t,x) clearly becomes delicate when f 1. In 
particular, direct use of the product (A.l) gives an algorithm that is only "linearly 
convergent" , i.e. the number of significant digits in the answer grows linearly with 
the number of terms taken. Moreover, the constant of proportionality in this relation 
is proportional to 1 — and thus deteriorates linearly as | 1. Finally, there is 
severe loss of numerical precision when multiplying numbers that are very near 1. An 
alternative approach can be based on the representation 

°° j-fc k 
k=l 

which is valid whenever \x\ < 1 and \tx\ < 1. This sum is again only linearly 
convergent, but the problem of loss of numerical precision is alleviated by use of the 
logarithm. 

A much more efficient algorithm can be based on the identity 

^ ('_t)na;n(n+l)/2 
R{t,x) = > ' „, r (A.3) 



n=0 
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due to Euler.^^ Because of the x"'^""'"^^/^ factor in the numerator, this algorithm is 
"quadratically convergent" : 



Proposition A.l ([61]) Define 



an — 



[1 - x){l - X^) ■ ■ ■ {1 - X"") 

Then, for \t\ < 1 and \x\ < e~'^ ('j > 0), we have 



(A.4) 



(a) An = 



(h) 5m = 



n=N 



< 



< 



g7r2/67-Ar(Ar+l)7/2 

1 _ e-(^+i)7 

g7rV37-A^(Af+l)7/2 



\R{t,x)\ - l-e-(^+i)^ 
Corolleiry A. 2 ([61]) Let K >0, and suppose that \t\ < 1 and \x\ < e~'^ ('j > 0). 



(a) IfN > 



7r2 2K 



+ , then Ajv < e 



-K 



37^ 7 



(b) IfN > J— + , then 5n < e"^ 



37^ 



7 



It turns out [61] that the a priori bound of Proposition A. 1(b) is asymptotically 
within 9.1% of being sharp when t = 1, < x = < 1 and N ^ l/'-f (moreover, 
in this case it is asymptotically sharp as 7 J, 0). But since this bound is overly 
pessimistic in other cases, it is of some value to provide an a posteriori bound on the 
truncation error that is more realistic, when x ^ (0, 1), than the a priori bound. Here 
is such a bound, which can be used a stopping criterion in the numerical algorithm: 

Proposition A.3 ([61]) Let \t\ < 1, \x\ < e"^ (7 > O; and N > (log2)/7. Then: 



(a) An = 



(b) S'n = 



n=N 
00 

n=N 



< \aN-i\ 



1 - 2e-^7 



< 



where Sn 



N-l 

E On 
n=0 



\Sm\ ~ \Sn\ l-2e-^T 
[Note also that 5n < 5^/(1 - 5^)-/ 
In particular, if N > (log3)/7; we have < |aiv-i| and S'j^ < |aAr-i|/|-S'jv|. 
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For a proof of (A.3), see e.g. [62, p. 19, Corollary 2.2], [63, p. 34, Lemma 4(a)] or [64, pp. 22-23]. 
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Let us conclude by making some brief remarks about the numerical precision that 
is required in intermediate stages of the calculation based on (A. 3). It turns out [61] 
that the largest term max |a„| can be as large in magnitude as e'^ /^^^ (and is indeed 

n 

of this order when < x < 1), while the answer R{t, x) can be as small in magnitude 
as g— ^ /67 jg indeed of this order when t — 1 and < a; < 1). It is there- 

fore necessary to maintain, in intermediate stages of the calculation, approximately 
(71^/47)/ log 10 ~ 1.07/7 digits of working precision beyond the number of significant 
digits desired in the final answer. 

We used all three algorithms — the product (A.l), the logarithmic sum (A. 2) and 
the quadratically convergent sum (A. 3) — and carefully cross-checked the value of 
R{t, x); we also verified numerically the error bounds of Proposition A.l, Corollary A. 2 
and Proposition A. 3. In order to guarantee that the roundoff error is under control, 
we performed all computations using Mathematica with a working precision of 
at least 100 digits and often much more (increasing the working precision until the 
answer is independent of the precision used) . 
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Lattice 


4th Zero 


5th Zero 


6th Zero 


3f 


X 3f 








3f 


X 6f 


2.552816126636 






3f 


X 9f 








3f 


X 12f 


2.562658027317 






3f 


X 15f 








3f 


X 18f 


2.565287184975 






3f 


X 21f 








3f 


X 24f 


2.566507072062 






3f 


X 27f 








3f 


X 30f 


2.567211365497 






4f 


X 4f 


2.604661945742 






4f 


X 8f 


2.618028652707 






4f 


X 12f 


2.618033986251 






4f 


X 16f 


2.618033988749 






4f 


X 20f 


2.618033988750 






4f 


X 24f 


2.618033988750 






4f 


X 28f 


2.618033988750 






4f 


X 32f 


2.618033988750 






4f 


X 36f 


2.618033988750 






4f 


X 40f 


2.618033988750 






5f 


X 5f 


2.618161303055 


2.795370504128 




5f 


X IOf 


2.618033988749 






5f 


X 15f 


2.618033988750 


2.947523648832 




5f 


X 20f 


2.618033988750 






5f 


X 25f 


2.618033988750 


2.968180058756 




5f 


X 30f 


2.618033988750 






5f 


X 35f 


2.618033988750 


2.976760450197 




5f 


X 40f 


2.618033988750 






5f 


X 45f 


2.618033988750 


2.981534673779 




5f 


X 50f 


2.618033988750 






6f 


X 6f 


2.618033979731 






6f 


X 12f 


2.618033988750 


3.001429148693 


3.054848659601 


6f 


X 18f 


2.618033988750 


3.000001523178 


3.100527321592 


6f 


X 24f 


2.618033988750 


3.000000001785 


3.118151997375 


6f 


X 30f 


2.618033988750 


3.000000000002 


3.127749140385 


6f 


X 36f 


2.618033988750 


3.000000000000 


3.133811079422 


6f 


X 42f 


2.618033988750 


3.000000000000 


3.137993327670 


6f 


X 48f 


2.618033988750 


3.000000000000 


3.141054810628 


6f 


X 54f 


z.oioOooyoo /^oO 


o. 000000000000 


0.1433900233 (a 


6f 


X 60f 


2.618033988750 


3.000000000000 


3.145239011028 


Beraha 


2.618033988750 


3 


3.246979603717 



Table 1: Real zeros of the chromatic polynomials of finite triangular-lattice strips with 
free boundary conditions in both directions, to 12 decimal places. A blank means that 
the zero in question is absent. The first three real zeros q = 0, 1, 2 are exact on all 
lattices. "Beraha" indicates the Beraha numbers E5 = (3 + a/5)/2, Bq = 3, and B^. 
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Lattice 


4tii Zero 


5tn Zero 


6tn Zero 


7tn Zero 


7f 


X 


7f 


2.618033988750 


2.978584823651 






7f 


X 


14f 


2.618033988750 


3.000000029690 


3.160410975706 




7f 


X 


21f 


2.618033988750 


3.000000000000 






7f 


X 


28f 


2.618033988750 


3.000000000000 


3.218685236695 




7f 


X 


35f 


2.618033988750 


3.000000000000 






7f 


X 


42f 


2.618033988750 


3.000000000000 


3.236121891966 




7f 


X 


49f 


2.618033988750 


3.000000000000 






7f 


X 


56f 


2.618033988750 


3.000000000000 


3.243833695579 




7f 


X 


63f 


2.618033988750 


3.000000000000 






7f 


X 


70f 


2.618033988750 


3.000000000000 


3.246633282347 




7f 


X 


77f 


2.618033988750 


3.000000000000 


3.247059872523 


3.254369173708 


7f 


X 


84f 


2.618033988750 


3.000000000000 


3.246965843358 




7f 


X 


91f 


2.618033988750 


3.000000000000 


3.246982133140 


3.258435734303 


8f 


X 


8f 


2.618033988750 


3.000359693703 


3.095706393163 




8f 


X 


16f 


2.618033988750 


3.000000000000 


3.229632685380 




8f 


X 


24f 


2.618033988750 


3.000000000000 


3.246928323759 




8f 


X 


32f 


2.618033988750 


3.000000000000 


3.246979586275 




8f 


X 


40f 


2.618033988750 


3.000000000000 


3.246979603712 




8f 


X 


48f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


56f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


64f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


72f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


80f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


88f 


2.618033988750 


3.000000000000 


3.246979603717 




8f 


X 


96f 


2.618033988750 


3.000000000000 


3.246979603717 




9f 


X 


9f 


2.618033988750 


2.999999518372 






9f 


X 


18f 


2.618033988750 


3.000000000000 


3.246969773686 




9f 


X 


27f 


2.618033988750 


3.000000000000 


3.246979603720 


3.342943823308 


9f 


X 


36f 


2.618033988750 


o r\r\r\r\r\r\r\r\r\r\r\r\ 

3.000000000000 


3.246979603717 




9f 


X 


45f 


z.OioUooyoo ( ou 


o.UUUUUUUUUUUU 




O.O / 4D4DZo4yo / 


9f 


X 


54f 


2.618033988750 


3.000000000000 


3.246979603717 




9f 


X 


63f 


2.618033988750 


3.000000000000 


3.246979603717 


3.387946181123 


9f 


X 


72f 


2.618033988750 


3.000000000000 


3.246979603717 




9f 


X 


81f 


2.618033988750 


3.000000000000 


3.246979603717 


3.395349738491 


9f 


X 


90f 


2.618033988750 


3.000000000000 


3.246979603717 




Beraha 


2.618033988750 


3 


3.246979603717 


3.414213562373 



Table 2: Real zeros of the chromatic polynomials of finite triangular- lattice strips 
with free boundary conditions in both directions, to 12 decimal places. We use the 
same notation as in Table 1. 
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Lattice 


4th Zero 


5th Zero 


6th Zero 


7th Zero 


8th Zero 


4p 


X 4p 


2.617986010522 


3 


3.465246100723 






A 

4p 


X op 


2.618033988740 


3 


3.475055224065 








X 12f 


2.618033988750 


3 


3.477452996799 






4p 


X IDp 


2.618033988750 


3 


3.478536268722 






4p 


X ZUp 


2.618033988750 


3 


3.479153472532 






4p 


X Z4p 


2.618033988750 


3 


3.479552148708 






/I 

4p 


X Zop 


2.618033988750 


3 


3.479830901859 






4p 


X OZp 


2.618033988750 


3 


3.480036768366 






4p 


X OOP 


2.618033988750 


3 


3.480195030232 






A ^ 
4p 


X 4Uf 


2.618033988750 


3 


3.480320488501 






5p 


X 5f 


2.618033990394 


3 








5p 


X lOp 


2.618033988750 


3 


3.196843987850 






Op 


X lop 


2.618033988750 


3 








Op 


X zUp 


2.618033988750 


3 


3.202699178454 






Op 


X zOp 


2.618033988750 


3 








op 


V Qn 
X oUp 


2.618033988750 


3 


3.204333275156 






F 

Op 


X oOp 


2.618033988750 


3 








Op 


X 4Uf 


2.618033988750 


3 


3.205100311429 






Op 


X 40f 


2.618033988750 


3 








op 


X OUp 


2.618033988750 


3 


3.205545558020 






6p 


X 6p 


2.618033988750 


3.001033705947 


3.125892136302 






6p 


X 12p 


2.618033988750 


3.000000003803 


3.198900652620 






6p 


X 18p 


2.618033988750 


3.000000000000 


3.217111179820 






Op 


X Z4p 


2.618033988750 


3.000000000000 


3.225649637432 






6p 


X 30p 


2.618033988750 


3.000000000000 


3.230657835149 






6p 


X 36p 


2.618033988750 


3.000000000000 


3.233968503481 






6p 


X 42p 


2.618033988750 


3.000000000000 


3.236327213212 






Dp 


X 4op 


2.618033988750 


3.000000000000 


3.238096251767 






Op 


X 04p 


2.618033988750 


3.000000000000 


3.239473538415 






Op 


X 60f 


2.618033988750 


3.000000000000 


3.240576619481 






Op 


X 66f 


2.618033988750 


3.000000000000 


3.241479828709 






Op 


X ( zp 


2.618033988750 


3.000000000000 


3.242232528364 






Dp 


X /op 


2.618033988750 


3.000000000000 


3.242868805497 






Dp 


X o4p 


2.618033988750 


3.000000000000 


3.243412961909 






Dp 


X yup 


2.618033988750 


3.000000000000 


3.243882786313 






7p 


X 7f 


2.618033988750 


3 


3.247001348628 


3.404690481534 




7p 


X 14f 


2.618033988750 


3 


3.246979603718 


3.414217072295 


3.458917430738 


'P 


X 21p 


2.618033988750 


3 


3.246979603717 


3.414213561735 




'P 


X Zop 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.470544903913 


' P 


X ODp 


2.618033988750 


3 


3.246979603717 


3.414213562373 




'P 


X 4zp 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.473634831556 


'P 


X 4yp 


2.618033988750 


3 


3.246979603717 


3.414213562373 




'P 


X OOp 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.475070205361 


'P 


X Oop 


2.618033988750 


3 


3.246979603717 


3.414213562373 




'P 


X (Up 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.475899672990 


8p 


X 8p 


2.618033988750 


3 


3.246979601854 






8p 


X 16p 


2.618033988750 


3 


3.246979603717 


3.414214415195 


3.472683999084 


8p 


X 24f 


2.618033988750 


3 


3.246979603717 


3.414213562387 


3.488644630018 


8p 


X 32p 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.495735217349 


8p 


X 40p 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.499773262291 


8p 


X 48f 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.502387969424 


8p 


X 56p 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.504221641913 


8p 


X 64f 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.505579831565 


8p 


X 72f 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.506626776159 


8p 


X SOp 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.507458740757 


Beraha 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.532088886238 



Table 3: Real zeros of the chromatic polynomials of finite triangular-lattice strips with 

periodic boundary conditions in the transverse direction and free boundary conditions 
in the longitudinal direction, to 12 decimal places. A blank means that the zero in 
question is absent. The first three real zeros q = 0,1,2 are exact on all lattices. 
"Beraha" indicates the Beraha numbers = {3 + \/E)/2, Bq = 3, Bf, Bg = 2 + \/2, 
and Bg. 
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Lattice 



4th Zero 



5th Zero 



6th Zero 



7th Zero 



8th Zero 



9th Zero 



10th Zero 



9p X 9f 
9p X 18p 
9p X 27p 
9p X 36p 
9p X 45p 
9p X 54p 
9p X 63f 
9p X 72f 
9p X 81p 
9p X 90f 



2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 



3.246980644227 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 



3.382733076359 

3.414215827400 
3.414213562359 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 



3.467483864312 
3.499429426359 
3.508825024982 
3.513393802382 

3.516109505154 



lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 
lOp X 



lOp 
20f 
30f 
40f 
50f 
60f 
70f 
80f 
90f 
lOOp 



2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 



3.414213601215 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 



3.522072913706 
3.532088885496 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 



3.618274945403 
3.618041117772 
3.618034257877 

3.618033998995 
3.618033989140 



3.620352727045 
3.623157797032 
3.624885448640 
3.626138707957 
3.627101535574 



lip 
lip 
lip 
lip 
lip 
lip 
lip 
lip 
lip 
lip 



lip 
22f 
33f 
44f 
55f 
66f 
77f 
88f 
99f 
IIOf 



2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 



3.414213539527 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 



3.532088885575 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 



3.608601511861 

3.617975980728 
3.618034451624 
3.618033985241 
3.618033988777 
3.618033988750 
3.618033988750 



3.627344614702 
3.631979435582 
3.634574709990 



12p X 
12p X 
12p X 

12p X 
12p X 
12p X 
12p X 
12p X 
12p X 
12p X 



12f 
24f 
36f 

48f 

60f 

72f 

84f 

96f 

108f 

120f 



2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 



3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 



3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 



3.414213593041 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 



3.511032635472 
3.532088885001 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 
3.532088886238 



3.618040035384 
3.618034035926 
3.618033989120 
3.618033988753 
3.618033988750 



3.624320958404 
3.627419917635 
3.629588218978 
3.631215401061 
3.632487726562 



Beraha 



2.618033988750 



3.246979603717 



3.414213562373 



3.532088886238 



3.618033988750 



3.682507065662 



Table 4: Real zeros of the chromatic polynomials of finite triangular-lattice strips with 
periodic boundary conditions in the transverse direction and free boundary conditions 
in the longitudinal direction, to 12 decimal places. We use the same notation as in 
Table 3. 
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Lattice 


4th Zero 


5th Zero 


6th Zero 


7th Zero 


8th Zero 


9th Zero 


10th Zero 


4z X 4f 
4z X 8f 

^Z ^ 

4z X 16p 
4z X 20p 
4z X 24p 
4z X 28p 
4z X 32p 
4z X 36p 
4z X 40p 


2.485072022789 
2.499965989337 
^.lyyyyyyo i ouo 
2.499999999885 
2.500000000000 
2.500000000000 
2.500000000000 
2.500000000000 
2.500000000000 
2.500000000000 


2.527537649962 
2.500034085574 

2.500000000115 

2.500000000000 
2.500000000000 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 


2.596617094656 
2.618031965217 

2.618033988750 

2.618033988750 
2.618033988750 










6z X 6f 
6z X 12p 

fir? Y 1 

6z X 24f 
6z X 30f 
6z X 36f 
6z X 42f 
6z X 48f 
6z X 54f 
6z X 60p 


2.618033988528 
2.618033988750 

2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 


3.000017186720 
3 0000000041 Q1 
3.000000000001 
3.000000000000 
2.722632835458 
2.722632835458 
2.722632835458 
2.722632835458 
2.722632835458 


3.117917986708 
3 1111 07899326 
3.150834657646 
3.156227017803 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 


3.159661924115 
3.162043675850 
3.163793040154 
3.165132700997 
3.166191662980 








8z X 8f 
8z X 16p 
8z X 24p 
8z X 32p 
8z X 40p 
8z X 48p 
8z X 56p 
8z X 64p 
8z X 72p 
8z X 80p 


2.618033988750 

2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 


3.000000844168 

3.000000000000 
3.000000000000 
2.821420495535 
2.821420495535 
2.821420495535 
2.821420495535 
2.821420495535 
2.821420495535 
2.821420495535 


3.203925019292 

3.246976356780 
3.246979603696 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 


3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 








lOz X IOf 
lOz X 20f 
lOz X 30f 
lOz X 40f 
lOz X 50f 
lOz X 60f 
lOz X 70f 
lOz X 80f 
lOz X 90f 
lOz X lOOp 


2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 
2.618033988750 


3.000000000000 
3.000000000000 
3.000000000000 
2.873731249334 
2.873731249334 
2.873731249334 
2.873731249334 
2.873731249334 
2.873731249334 
2.873731249334 


3.246953122227 
3.246979603717 
3.246979603717 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 
3.000000000000 


3.416320582746 
3.414213564771 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 
3.246979603717 


3.428202969384 
3.470075808656 
3.414213562373 
3.414213562373 
3.414213562373 
3.383128531235 
3.383128531235 
3.383128531235 
3.383128531235 


3.484585415185 
3.492475877808 
3.497477600415 
3.414213562373 
3.414213562373 
3.414213562373 
3.414213562373 


3.500942791087 
3.503489509874 
3.505442176204 
3.506987965042 


Beraha 


2.618033988750 


3 


3.246979603717 


3.414213562373 


3.532088886238 


3.618033988750 


3.682507065662 



Table 5: Real zeros of the chromatic polynomials of finite triangular-lattice strips 
with zig-zag boundary conditions, to 12 decimal places. We use the same notation as 
in Table 3. 
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Eigouvahu 


^-Crossing C 






Isolaic'd Points 


Lattice 


#C 


#E 


#T 


# D 


# ER 


min Re q 




max Re q 


#RI 


# CI 


2f 


















2 





3f 


3 


6 











1.20474 


2.56984 


3.40223 


3 





4f 

r 


2 


6 


2 








0.81647 


2.75925 it 0.15444 z* 


3.63983 


4 





5f 


1 


6 


12 





4 


0.55862 


3 


3.77830 


4 





6 k 


1 


6 


4 








0.37963 


3.16093 


3.86641 


5 





7f 


it 


6t 


4t 


n 


ot 


0.25054 


3.27640 


3.92580 


R 


n 


8f 


it 


6t 


4t 


n 


ot 


0.13343 


3.36106 


3.96756 


a 


n 


9f 














3.42513 




7 





Zp 


















2 


n 


3p 


















3 





4p 


3 


4 








1 


1.37053 


3.48141 


4 


5 





5p 


3 


4 





n 


1 

± 


0.47725 


3.20722 


3.87699 


r. 


n 


6p 


1 


4 


2 








0.02077 


3.25242 


4.28386 


6 





7p 


1 


4 


4 





1 


-0.22792 


3.47900 


3.99964 


7 





8p 


1 


4 


2 








-0.37137 


3.51477 


4.04970 


7 





9p 


1 


4 


2 








-0.45760 


3.52706 


4.28286 


7 





lOp 


It 


4t 


2t 





ot 


-0.51081 


3.63483 


4.12341 


9 





lip 


It 


4t 


2t 





ot 


-0.54399 


3.64414 


4.15609 


9 





12p 






2t 








3.64317 




9 





2z 


















2 





4z 


2 


4 











2.09914 


2.73717 ± 0.17233 i* 


4.00485 


5 





6z 


1 


4 


2 








0.36185 


3.17526 


4.25895 


6 


1 


8z 


It 


4t 


2t 





ot 


-0.21435 


3.39410 


4.28991 


7 


It 


lOz 














3.52044 




9 





Table 6: Summary of qualitative results for the eigenvalue-crossing curves B and for 
the isolated limiting points of zeros. For each triangular-lattice strip considered in 
this paper, we give the number of connected components of B (7^ C), the number 
of endpoints (7^ E), the number of T points (7^ T), the number of double points (7^ 
D), and the number of enclosed regions (7^ ER-); we also give the minimum value 
of Req on B, the smallest value qo where B intersects the real axis (* denotes an 
almost-crossing), and the maximum value of Reg on B. We also report the number 
of real isolated limiting points of zeros Rl) and the number of complex-conjugate 
pairs of isolated limiting points (# CI). The symbol t indicates uncertain results. 



64 



m 


TriCyl(m) 


# VA 


TriCyl'(m) 


SqCyl(m) 


1 


1 





1 


1 


2 


1 





1 


1 


3 


1 





1 


1 


4 


2 





2 


2 


5 


2 





2 


2 


6 


5 





5 


5 


7 


6 





6 


6 


8 


15 


2 


13 


14 


9 


28 


12 


16 


22 


10 


67 


32 


35 


51 


11 


145 


100 


45 


95 


12 


368 


272 


96 


232 


13 


870 


742* 


126* 


498 


14 


2211 


1940* 


267* 


1239 



Table 7: Transfer-matrix dimensions for a triangular-lattice strip of width m and 
cylindrical boundary conditions. For each value of the strip width m we give the 
dimension of the transfer matrix [TriCyl(m)], the number of vanishing amplitudes 
VA), and the effective dimension of the transfer matrix [TriCyl'(m) = TriCyl(m) — 
# VA]. For comparison, we also give the dimensionahty of the transfer matrix for a 
square-lattice strip of width m and cylindrical boundary conditions [SqCyl(m)]. The 
values of TriCyl(m) and SqCyl(m) were obtained in [16, 53]. An asterisk denotes 
conjectured results. 
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Table 8: Endpoints in the complex p-plane with |p| = 1 of the curves where \g2/gi\ = 
1. For each endpoint we show the "Curve" to which it belongs (see Figure 22), the 
estimated value of its phase Argppit (see text), the conjectured exact value Argp, the 
"discrepancy" [cf. (6.10)], the "significance level" e [cf. (6.14)], and the corresponding 
"primary" 9 value. For some selected values we also include a label ("Point"). When 
we have a curve for which one endpoint is well-determined and the other is not, we 
include the former point in the category "Others" (see text). 
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Table 9: Endpoints in the complex y- plane with |y| = 1 of the curves where \gz/g\\ — 
1. For each endpoint we show the "Curve" to which it belongs (see Figure 23), the 
estimated value of its phase Argypit (see text), the conjectured exact value Argp, the 
"discrepancy" [cf. (6.10)], the "significance level" e [cf. (6.14)], and the corresponding 
"primary" 6 value. For some selected values we also include a label ("Point"). When 
we have a curve for which one endpoint is well-determined and the other is not, we 
include the former point in the category "Others" (see text) . 
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Table 10: Beraha numbers Bn that are zeros of det D{q). Those shown in boldface 
(resp. normal face) correspond to the vanishing of a dominant (resp. subdominant) 
amplitude. 
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(a) 




(b) 

Figure 1: Two ways of building a triangular- lattice strip using a transfer- matrix 
approach, (a) Standard method (see e.g. Ref. [16]). (b) Alternative method (called 
"zig-zag" boundary conditions). 
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Figure 2: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices Sp x 15^? (squares), Sp x 30^? (circles) and 3^ x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 3: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices x 20^? (squares), x AOp (circles) and x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 4: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices Si;' x 25^? (squares), 5f x 50p (circles) and dp x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 5: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices Gp x 30^? (squares), x 60^? (circles) and x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 6: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 7f x 35p (squares), 7f x 70p (circles) and x oop (solid line). 
The isolated limiting points are denoted by a x . 



74 




2 4 

Re(q) 

Figure 7: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices Si;' x 40^? (squares), Sp x 80p (circles) and x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 8: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices dp x 45^? (squares) and 9p x 90 p (circles). The isolated limiting 
points are denoted by a x. 



76 



Zeros tri lattice = 4p 

' \ ' r 




I I I I I I I 

2 4 

Re(q) 

Figure 9: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 4p x 20^? (squares), 4p x AOp (circles) and 4p x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 10: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 5p x 25^? (squares), 5p x 50^ (circles) and 5p x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 11: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 6p x 30^? (squares), 6p x QQp (circles) and Qp x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 12: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 7p x 35^? (squares), 7p x 70p (circles) and 7p x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 13: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 8p x 40^? (squares), Bp x 80p (circles) and 8p x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 14: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 9p x 45^? (squares), 9p x 90^ (circles) and 9p x oop (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 15: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices lOpxbOp (squares), lOp x lOOp (circles) and lOp x cxd^ (sohd line). 
The isolated limiting points are denoted by a x . 
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Figure 16: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices lip x SSi? (squares), lip x llOi;' (circles) and lip x cxDp (sohd line). 
The isolated limiting points are denoted by a x . 
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Figure 17: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 12p x GOi? (squares) and 12p x 120^ (circles). The isolated limiting 
points are denoted by a x. 
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Figure 18: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 4p x 20 z (squares), 4p x 40 z (circles) and 4p x ooz (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 19: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 6p x 30 z (squares), 6p x 60 z (circles) and 6p x ooz (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 20: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices 8p x 40 z (squares), Bp x 80 z (circles) and 8p x ooz (solid line). 
The isolated limiting points are denoted by a x . 
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Figure 21: Zeros of the partition function of the g-state Potts antiferromagnet on the 
triangular lattices lOp x 50 z (squares) and lOp x 80^ (circles). The isolated limiting 
points are denoted by a x. 
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Igg/gJ = 1 in complex p— plane 
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Figure 22: Equimodular curves |5'2/5'i| = 1 in the complex p-plane. Several important 
points are labelled G, H, I, J, . . . (see text). 
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Figure 23: Equimodular curves \gzl g\\ = 1 in the complex y-plane. Several important 
points are labelled G, H, I, J, . . . (see text). 
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Figure 24: Equimodular curves for the eigenvalues gi in the complex ^-plane. The 
portions of curves where the equimodular eigenvalues are dominant (resp. subdomi- 
nant) are depicted in black (resp. pink). The eigenvalue gi is dominant in each region 
labelled Pj. Several important points are labelled A, B, C, . . . (see text). To facili- 
tate comparison with Baxter's results [19, Fig. 5], we have used the same labelling of 
points wherever possible. 
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Figure 25: Dominant equimodular curves for the eigenvalues gi in the complex q- 
plane (in black). For comparison, we show (in pink) the limiting curve B for the strip 
= lip. To facilitate comparison with Baxter's results [19, Fig. 5], we have used 
the same labelling of points wherever possible. We warn the reader that the presence 
of the additional curves (AH, LK, etc.) to the right of g = has not been definitively 
established; see Sections 6.3 and 6.4 for a detailed discussion. 
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Limiting Curves Tri Lattice 




3.6 3.8 4 4.2 

Re(q) 

Figure 26: Detail of Figure 25 around the point g = 4. We depict the dominant 
(resp. subdominant) equimodular curves in black (resp. pink). The solid circles (•) 
denote special points discussed in the text, while the squares (□) and empty circles 
(o) denote the Beraha numbers q = Bn, . . . ,Biq, B24, -B27 and B^o = 4 (point C). 
The empty circles denote those Beraha numbers which belong to any equimodular 
curve. We denote by the regions where the eigenvalue Qi is dominant. 
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Figure 27: Limiting curves for the triangular-lattice strips Lp x oop with 3 < L < 8. 
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Figure 28: Limiting curves for the triangular-lattice strips Lp x cxDp with 4 < L < 11. 
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Figure 29: Limiting curves for the triangular-lattice strips Li x oof with L = 4,6, 8. 
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AT tri-lattice Potts Im(q)=2 




Re(q) 

Figure 30: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 2. Curves show the logarithm of Baxter's 
eigenvalues log \gi\: the solid curve corresponds to gi, the dashed curve to g2, and the 
dot-dashed curve to g^. Points show the free energy fk = (l/L)log|Afc| associated 
to the five largest eigenvalues (in modulus) of the transfer matrix for a triangular- 
lattice strip of width L = 12p. The index k is coded as follows: k = 1 (the dominant 
eigenvalue, o), 2 (□), 3 (A), 4 (o), and 5 (x). The solid green dots represent the 
extrapolation of the dominant eigenvalue to L = oo. 
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AF tri-lattice Potts Im(q)=1.5 
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Figure 31: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 1.5. The notation is as in Figure 30. 
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AF tri-lattice Potts Im(q)=l 
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Figure 32: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 1. The notation is as in Figure 30. 
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AF tri-lattice Potts Im(q)=0.5 
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Figure 33: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 0.5. The notation is as in Figure 30. 



101 



AF tri-lattice Potts Im(q)=0.35 
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Figure 34: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 0.35. The notation is as in Figure 30. 
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AF tri-lattice Potts Im(q)=0.2 
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Figure 35: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 0.2. The notation is as in Figure 30. 
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AF tri-lattice Potts Im(q)=0.1 
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Figure 36: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 0.1. The notation is as in Figure 30. 
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AF tri-lattice Potts lm(q)=0 
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Figure 37: Comparison of Baxter's eigenvalues (6.2) with the most dominant eigen- 
values of the transfer matrix for Img = 0. The notation is as in Figure 30. Baxter's 
eigenvalues g2 [cf. (6.2b)] and g-i [cf. (6.2c)] are depicted only on the interval < g < 4, 
while gi [cf. (6.2a)] is shown only outside that interval (see text for details). 
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